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Series Editor's Introduction 


R. Robert Huckfeldt, C. W. Kohfeld, and Thomas W. Likens have 
added a new dimension to this series by writing Dynamic Modeling: An 
Introduction. In this monograph, they introduce the reader to the tools 
necessary for sophisticated modeling of the processes of social change, 
broadly understood. Their emphasis is on the use of difference equations in 
modeling such changes, although their treatment of the topic is such that 
no advanced training in mathematics is required. Unlike differential 
equation models, in which time is treated as continuous and rates of change 
are considered to be instantaneous, in difference equation models, time is 
treated asa series of discrete, equally spaced units. The latter have probably 
been given less attention by social scientists, but in many instances, such 
models often make more substantive sense. 

Using examples that include models of social and political mobilization, 
diffusion, and social interaction, Huckfeldt et al. introduce the reader to 
first-order linear and nonlinear difference equation models in Chapters 
2 and 3. They show the novice how to examine the conditions of equi- 
librium and stability, and how to estimate the parameters of each model. 
Examples abound, simplifying the presentation and assisting the mathe- 
matically untrained reader in his or her attempt to understand the appli- 
cation of difference equation models in the social sciences. 

In Chapters 4 and 5, the authors move beyond first-order equations, 
discussing linear and nonlinear higher-order equations in the two chapters, 
respectively. There is an extended example of arms race models, developed 
by Richardson and expanded upon by Gillespie and Zinnes in recent 
years, as well as modeling examples of legislative review and of the 
budgetary process. 

With exception of Chapter 5, the authors have presented this reason- 
ably complex material without assuming much more than high school 
algebra. The mathematics developed in Chapters 5, on nonlinear dynamic 
systems, is more demanding but the authors’ intention “is not to equip 
readers fully with the analytical methods required to conduct similar 
analyses, but rather to illustrate the direction in which the work in this 
monograph naturally moves.” The extended example of budgetary com- 
petition is therefore meant as illustrative and should not be beyond the 
reach of most of our readers when they read it within the framework of 
the authors’ intentions. 

Most of the previous monographs published in this series have been 
purely statistical in orientation. This monograph is an introduction to 
mathematical modeling and as such represents the first in a series of 
departures from that theme. The Quantitative Applications series will 
shortly publish monographs on nonstatistical topics such as questionnaire 
design and construction. As the topics covered by this series broaden, our 
intentions remain firm: to present simple, readable introductions to topics 
of interest to social scientists who engage primarily in quantitative work. 


—John L. Sullivan 
Series Co-Editor 
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1. DYNAMIC MODELS AND SOCIAL CHANGE 


Social reality is more than a set of substantively related but randomly 
occurring events. In this monograph we think of social phenomena as 
processes—structured series of events, operations, and activities whose 
logic is orderly and predictable. The goal of dynamic modeling is to 
specify the structure of such processes and to deduce the manner in which 
they generate social change. 

This monograph presents a set of techniques and strategies for dynamic 
modeling within the context of several substantive applications. One 
component demonstrates the translation of substantive ideas into precise, 
mathematical statements. We are concerned with a technology for speci- 
fying models that have (1) mathematically based substantive interpre- 
tations, and (2) logically tractable and meaningful deductive consequences. 
Indeed, descriptive adequacy—a persuasive descriptive interpretation for 
a model and its consequences—is the main technical criterion for model 
evaluation. 


Difference Equations as Representations of Change 


The notion of change through time iscentralto the modeling techniques 
developed here, and the reader is familiarized with the use of difference 
equations to represent processes of change. Difference equation models 
rest upon well-developed mathematical theory, enabling a straight- 
forward and complete analysis of a process through time. They do not, 
however, require advanced mathematical training to be understood. A 
large majority of what follows requires only that readers possess a working 
knowledge of high school algebra and, perhaps, some patience. 


Difference equation models are not the only dynamic models that can 
represent change in some phenomenon over time. Time can be treated as 
continuous, thereby producing differential equation models based on 
infinitely small time units and instantaneous rates of change. In difference 
equations time is treated as a series of discrete, equally spaced units. The 
choice between analogous pairs of difference and differential equation 
models is easily made when phenomena have a natural time metric in 
either discrete or continuous time. There are often grounds, however, for 
viewing discrete events as manifestations of an underlying continuous 
process, or continuous processes as being manifested at discrete points in 
time, and hence the argument that phenomena have a "natural" repre- 
sentation in discrete or continuous time is frequently undermined. Indeed, 
arguments can be made to justify representing most processes as either 
discrete or continuous (Kohfeld and Salert, 1982). 

Several considerations lead us to develop dynamic models as discrete 
time representations. First, observations always occur in discrete time. 
In the physical sciences observations can be taken in such rapid succession 
that they approach a continuous observation, but this is seldom the case 
in the social sciences. Data on social phenomena are typically available 

at discrete time points that are significantly separated in time. For example, 
a scientist who is interested in the cooling process of a heated liquid takes 
thousands of measurements separated only by small fractions of a second. 
A social scientist who is interested in the socialization process of a young 
child acquires far fewer measurements that are much farther apart in time 
As a result, we frequently base our representation on the observed behavior 
rather than the underlying process. 

Second, given analogous difference and differential equations, differ- 
ence equation models are capable of representing a wider vatiet f 
possible behavior patterns over time, and thus offer greater Фан ol 
for modeling observed behaviors with simpler models. This feat R 
difference equations is convenient, but also mildly disturbing. | ature of 
that the choice of a modeling strategy and an accompanying r В. USUBgests 
for time is not merely a technical matter unrelated to р entation 
siderations. Rather, the choice of discrete time modes s stantive con- 
tinuous time models has important substantive im See! than con- 

Appendix for a more complete discussion of these te ications (see the 

Finally, and not insignificantly, difference equati atters), 
lower level of mathematical sophistication, We co s models require a 
the same ground in the available space for the dme uld not hope to cover 
ment proceeded in terms of differential equation readers if our de 

We embrace the notion that any model pia se 
representation of a process. Even in representin ae only as an abstract 
employ an abstraction with inevitable limitatio B time we are forced to 
here should not be judged on the extent to are The Models employed 
process, but rather on the extent to which they "e they replicate a social 
process. 


velop- 


Pas Understand a social 


Synchronic versus Diachronic Change 


Throughout this monograph we utilize the concept of the structure of 
a dynamic process. By this we mean the laws of change that define the 
time-dependence of a process. The first objective of dynamic modeling is 
to understand, substantively, the mechanisms that are generating change 
in some observable phenomenon, and then to translate this set of ideas 
about structure into a mathematical language that can be deductively 
investigated. 

The techniques presented here are intended to equip the reader with a 
set of tools for representing dynamic structures and investigating how such 
structures behave. It must be emphasized, however, that these technologies 
are only applicable to synchronic change—to change produced by a fixed 
dynamic structure. The treatment of diachronic change—the process of 
moving from one structure into another—requires different method- 
ologies and is not considered here. 

In Chapter 2, for example, we present a model of mobilization. We 
hypothesize a very elementary dynamic structure by which mobilization 
levels change from one point in time to the next. We use a difference 
equation model to represent this structure, and to study its properties. 
The problem of how one mobilization process is transformed into another, 
however, cannot be addressed within this technology. Our modeling allows 
us to learn essentially all there is to know about any particular mobilization 
process; it may even help to understand the conditions that will cause 
one process to break down and be replaced by another. But it will not 
allow us to deduce what new structure will emerge. Theories attempting 
to explain the origin of a new dynamic structure, or the replacement of 
One structure by another, are beyond the bounds of the models considered 
here? 


Dynamic Models and Other Models of Change 


Recent developments in social science research have more and more 
come to include time as an important explanatory concept. Dynamic 
modeling is only one means for addressing the importance of time and 
change, and it is valuable to juxtapose dynamic models with other models 
and approaches for understanding social processes and social change. 

Various statistical strategies, and particularly statistical time series 
analysis, are frequently used to incorporate time within social science 
explanations. For example, Tufte (1978) uses a statistical strategy to 
relate economic conditions in the spring to election outcomes in the 
subsequent fall. Lewis-Beck and Alford (1980) use a time series design to 
relate policy interventions in coal mines with changes over time in mine 
safety. And, Hibbs (1976) uses a different time series design to show the 
over-time relationship between unemployment rates and labor union 
militance. Many of these efforts rest upon a well-developed econometrics 


10 


literature that is concerned with a variety of estimation problems. Indeed. 
a continuing line of inquiry is related to many of the estimation issues 
that arise within time series designs (Ostrom, 1978). 

Dynamic models are complementary to these statistical approaches, but 
they are distinct from them as well. This is not to say that issues in statistical 
estimation are irrelevant to dynamic models: One strategy for obtaining 
values for structural parameters within a dynamic model is to estimate 
them statistically. Indeed, during the course of our development we 
statistically estimate model parameters in several different instances. The 
important point remains: Dynamic modeling is not a statistical technique. 
and this monograph is not statistical in content. 


The Structure of the Presentation 


This monograph considers four general types of increasingly complex 
dynamic models in the ensuing four chapters. These chapters deal. 
respectively, with single-equation linear models, single-equation nonlinear 
models, linear systems composed of linear interdependent equations, and 
nonlinear systems composed of nonlinear interdependent equations. The 
substantive topics of these models— political mobilization, social dif- 
fusion, arms races, and competition over scarce resources--are funda- 
mental to a variety of social science concerns but incidental to the objective 
of the manuscript. The substantive topics are used essentially as vehicles 
for introducing a series of issues related to dynamic modeling in discrete 
time, and these issues are germane to an even broader range of substantive 
concerns. 

In summary, our goal is to present a framework for considering the 
logical consequences of a social process with a particular and specified 
structure. The reader who perseveres will acquire an overview of four 
broad classes of models and an understanding of the research situations 
in which they are substantively appropriate and likely to be useful In 
addition, we hope that some specific modeling Strategies will be ac wired 
in the form of mathematical representations for frequently recurrin x ial 
phenomena. We begin our task by developing a dynamic wad ртт 
political mobilization process. el of the 


2. FIRST-ORDER LINEAR DIFFE 
MODELS: A MODEL OF THE MOBILIZ 4 ЕОЧАТОМ 


ILIZATION PROCESS 
This chapter examines some element i 

mobilization and develops a dynamic A ca re 
We begin with a set of statements about mobilizati mo 
statements are then translated into a formal, e a mechanisms, These 
Finally, techniques for analyzing and estimating th matical representation. 
During the course of developing this mad 8 the model are introduced. 


el, wi i 
we define ang illustrate a 


garding political 
bilization process. 
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number of ideas that are central to dynamic modeling. These include 
difference equations and their solutions, the equilibrium of a dynamic 
process, system stability, and qualitative behaviors over time. By intro- 
ducing these concepts within the context of a particular model, the power 
of dynamic technologies as a strategy for investigating social and political 
problems becomes apparent. 


Mobilization Processes 


Social and political events occur at specific times and in specific con- 
texts, but they are seldom divorced from the past or from the future. Many 
of these events are imbedded within processes that can be best understood 
as occurring across time. Political or social mobilization—the recruitment 
of individuals to some cause (or party or belief or practice)—is one such 
phenomenon. A current mobilization level is inherently dependent upon 
mobilization in the past, just as future mobilization is tied to mobilization 
in the present. In short, the process of mobilization is time-dependent in 
the sense of being sequenced. Future levels of mobilization are intimately 
connected with past levels according to a systematic dynamic structure. 
Using a dynamic model, then, we wish to develop a specific representation 
of the logic underlying change in mobilization levels from one point in 
time to the next. 


A First Dynamic Model: The Gain/Loss Formulation 


Any population can be usefully partitioned into two groups relative to 
the mobilization process: those individuals who are currently mobilized 
and those who are not. Denote as M, the proportion of the relevant popu- 
lation which is mobilized at time t. Since the partitioning is exhaustive, the 
proportion of the population that is not mobilized at time t is given by the 
quantity | - Mi. М, is in principle observable, and it is formally known as 
a system state—a numeric quantity characterizing a process at a point 
in time. 

The problem is to model the change in M, across discrete points in time. 
Change in M, may be written formally as: 


AM: = Mer - Mı [2.1] 


where 7 is the differencing interval. For our treatment of discrete-time 
dynamics, the differencing interval is always assumed to be unity unless 
otherwise noted. While the assumption of 7 = l is a standard modeling 
practice, the essential point for technical development is that 7 is fixed. 
A linear transformation can always be used to obtain a 7* that has the 
value unity and similarly, the inverse transformation can be used to recover 
the interval in its original units. We also follow the usual practice of 
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developing difference equation dynamics as functions on the integers. 


but the key formal properties are a fixed difference interval and a well- 
defined, uninterrupted sequence. 


We thus define the first difference of M, as: 


AM: = Ma - М, pa 


Higher-order differences obtained by continued differencing of the 
dynamic equation are considered in later chapters. 

A change in the level of political mobilization, AM,, results from two 
simultaneous sources. From time t to (t+ 1) some individuals who are not 
mobilized will be recruited, generating a gain in Mı. Some other individuals 
who have already been mobilized will defect, thus producing a /oss in M. 
The net effect of these simultaneous (but independent) gains and losses 
gives the change in М,, AM.. 

Let a constant parameter, 8, represent the probability that an un- 
mobilized individual is mobilized from time t to (t + 1). Then the gain in 


М, is given by the recruitment rate, в, multiplied times the unmobilized 
population, (1 - M): 


gain of recruits = g(1 - М) [2.3] 


Similarly, define a second constant parameter, f, to represent the prob- 
ability that an individual who is mobilized at t will case to be mobilized 
at(t+ 1), The loss in M, is calculated as this defection rate multiplied times 
the mobilized population: 


loss of members - fM, [2.4] 


Thus, change in mobilization across time is given by: 


AM. = g(1 - M) - fM, [2.5] 


(Change = gains – losses) 


Equation 2.5 is the sim 


i í of gain/loss models 
whose underlying logic expresses the change in some stat : tem as 
the sum of gains and losses that have Occurred € of a syste 


changes his or her behavior. This ilisti 

r. Probabilistic ; ati 
course, change the deterministic nature of the aioe oat ca 
whether one thinks of g and f as rates oras Probabili i "i Мше 
are logically restricted to the (0,1) interval 54 noe 
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Within these (0,1) bounds, the magnitudes of both g and f depend upon 
the time metric chosen to represent a process. This dependence is easily 
illustrated by comparing presidential popularity to peasant farmers’ 
adoption of agricultural innovations. Investigators may focus on short- 
term fluctuations (7 of one month), or longer-term fluctuations (7 of one 
year). Months might be more appropriate for presidential popularity, and 
years might be more appropriate for changes in peasant agricultural 
practices. 

The descriptive adequacy of equation 2.5 is considerably improved by 
the incorporation of one additional idea. For many mobilization processes, 
some individuals are simply not susceptible to recruitment efforts under 
any circumstances. Our model accommodates this notion by letting an 
additional parameter, W, represent the proportion of individuals in the 
population who cannot or will not be mobilized. We wish to remove those 
individuals from consideration, and so the term | - W gives the maximum 
mobilization level possible. At any particular moment, then, the popula- 
tion susceptible to recruitment is given by the quantity | - W - M.. It is 
convenient, however, to define the quantity. 


рес [2.6] 


such that L may be interpreted as the upper limit on the potentially 
mobilized population. Substituting, our model becomes: 


AM, = g(L - М) - fM, [2.7] 


The parameter W (ог 1.) has a substantive and logical status somewhat 
different from that of g and f—the gain and loss rate parameters. The 
quantities g and f are rates or probabilities, while the parameter W (or L) 
is measured in the metric of the state M, and constitutes a modification 
in the description of the population that is subject to the mobilization 
process. In a purely formal sense, g, f, and W are all parameters of the 
process, but in a substantive sense, g and f have a different epistemological 
status when contrasted with W. Put another way, W describes while g and 
f give laws of change. 

The model treats the population of interest as three distinct groups: 
the mobilized (M,), the unmobilized subject to recruitment (L - M;), and 
those who cannot be mobilized (W = | - L). The state variable M, is ob- 
servable and subject to variation across time. The parameters f, g, and L 
are assumed to be invariant over time; they are лог directly observable 
and must be estimated empirically. 

Throughout this chapter the gain/loss model is discussed in terms of 
political parties and their electoral support.This convention facilitates 
our discussion, but should not obscure the broad substantive applicability 
of the model. This particular model, and others similar to it, are useful in 
understanding a wide range of social and political phenomena: support 
for political candidates during a campaign, changes in consumer prefer- 
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ence, adoption of birth control practices, disapproval of the Vietnam War. 
and so forth. 


Difference Equations as Representations of Change 

Equation 2.7 expresses the /ogic of the model most directly. In analyzing 
difference equation models, however, it is convenient to rearrange the 
original substantive specification algebraically in order to produce the 
following standard form: 


Mei = ao + aiMi [2.8] 


This is readily accomplished by recalling that 


AM: = Mui - M [2.9] 
Substituting equation 2.9 into equation 2.7 and expanding yields 

Mui = Mi = gL- gM: - fM, [2.10] 
or, regrouping and collecting terms: 

Mui = gh + (1 - f- gM [2.11] 


Тһе геайег should verify that equation 2.11 is in the standard form for 
first-order linear difference equations given by equation 2.8, with 


ao = gL [2.12] 


[2.13] 


a-l-f-g 


The standard form of equation 2.11 is crucialt 
mobilization model's mathematical and deductive properties, W. 
terize this form as a first-order linear difference Satis kx ; € charac- 
coefficients. The dynamics are expressed as a difference RE wit п constant 
time is taken as a series of discrete, equally spaced mE Бера 
system states are sequenced on time (i.e., they are time a and because 
equation is /inear, because none of the a, Coefficients ar ependent). The 
of the states of the system (the Mi, fork = 0,1,2 ) NS function of any 
tinction between "nonlinear" and "linear" difference 6 їп bn hat tepis: 
to the shape of the generated time path for the EU isunrelated 
typically produce nonlinear paths, as do nonlinear ig Linear models 

The equation has constant coefficients because th сов 
fixed model parameters (g, f, L) that do not ry RON are functions of 
simplifying assumption. Finally, the order of а duel time—a typical 
the maximum number of lags required to generate тепсе equation gives 


any particular Observa- 


oan understanding of the 


is 


tion in time; it is determined by subtracting the lowest time subscript in 
the equation from the highest time subscript. In this instance only one 
time point in the past is required to predict the next point, and hence the 
model is first order: 


(pelystst [2.14] 


More complex difference equations are addressed in subsequent chapters. 
As we shall see, however, the first-order linear difference equation with 
constant coefficients is a powerful tool for analyzing and describing laws 
of change. 

The next task is to ascertain the specific deductive consequences of the 
intuitive ideas regarding mobilization that are expressed in the model. 
What does the model assert, substantively, about mobilization processes? 
What histories of mobilization does it predict, under what circumstances? 
How well do these predictions stand against data? In order to address 
these issues, we briefly introduce the idea of so/utions to difference equa- 
tions and the techniques by which such solutions are derived. 


Solutions to First-Order Difference Equations 


A solution to a difference equation is a single function that generates a 
sequence of values satisfying the difference equation at each time point 
(Goldberg, 1958). The solution for equation 2.8 is developed in a straight- 
forward fashion by first generating a series of equations that calculate 
several successive states as functions of the initial condition. 


ам 


*aM,- ag ajay tay Mo 


2 3 
ag tay dy t ajag t a M 


My =a *24M, , 
2. n-l, Qn 
Sag taag арар t sss жар ap +аМу 
n 
: [2.15] 
Mn = ao [2.16] 


for all time. Thus, the process is stationary, and the future is constant. 
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Alternatively, suppose that ao = 0 and invoke equation 2.15 once again. 
In this instance the process is time dependent and defined as 
M, = a]Mo [2.17] 
These two hypothesized conditions (ao = 0 and a; = 0) demonstrate that a 
generated sequence will fail to be dynamically interesting unless ay 15 
different from zero. (Although, as we will see below, it is still possible for 
ài to be different from zero and for the generated sequence to be constant 
over time.) 

Equation 2.15 fulfills the definition for a solution to a difference equa- 
tion by generating a sequence that satisfies the difference equation at each 
time point, but it offers little advantage over the alternative method for 
calculating M,—successive calculations using the original dynamic 
equation. More importantly, it is not analytically useful. The solution 


becomes more instructive when it is translated into its closed form. 
Begin by expressing equation 2.15 as 


M. =S + a]Ms [2.18] 
where 

5 = ao + aiao + afao +... + alla, [2.19] 
Multiply both sides of equation 2.19 by a, to produce 

aS = aiao + afao +... + atag [2.20] 


Finally, subtract equation 2.20 fr. 


) от equation 2.19, and divide both sides 
by | - ai, given that ài * 1, to obtain 


This solution is clearl 


y degenerate when a, 
of zero in the fi 


| = 1, producing a denominator 
116 list term on the right-hand side. In this special case (ai = 1) 
the following closed-form solution is readily derived: 

Ms = ao + aiMs-i = nay + My [2.23] 
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Summarizing our results, first-order linear difference equations with 
constant coefficients, written in the standard form (from equation 2.8): 


M, = ao + а, Ма. 


have solutions given by (from equations 2.22 and 2.23, respectively): 


= 
D] 


+alMy (fora, #1) 


= 
[ 


= nay, tM, (fora = 1) 


Furthermore, it can be shown (Goldberg, 1958) that these are the only 
solutions. Thus, given any dynamic model written in the form of a first- 
order linear difference equation with constant coefficients, we are able 
to obtain a solution; and once the parameters ao and a; are known, we can 
directly ascertain the value of the time series at any n" point in time. 

The utility of the solution goes far beyond merely calculating sequence 
values. We wish to interpret the solution substantively, and it is to this 
substantive potential that we now turn. 


Equilibrium and Stability 


Before treating equilibrium in its technical aspects, several substantive 
examples are in order. The general level of partisan support for a political 
party in the United States exhibits remarkable stability over long time 
periods, often referred to as political epochs (Burnham, 1970). These 
periods are typically ushered in by an election or series of elections in 
which a rapid change in relative party fortunes is registered. This well- 
known phenomenon in American electoral history is an example of a 
dynamic system that experiences a fundamental alteration: an important 
parameter or input changes, producing rapid adjustment of the process 
to a new equilibrium. Equilibrium does not mean the absence of change in 
individual behavior but rather the typical steadiness of a measure of 
behavior. That is to say, equilibrium as a technical concept refers to a 
numerical measure of the level of the process. Clearly there are an un- 
limited number of ways in which the Democratic party might receive 58% 
of the presidential vote from election to election. People die, new electors 
join, citizenship is granted to immigrants, general levels of partisanship 
may change, yet the percentage may appear nearly constant. It is likely 
that this is a substantive example of equilibrium. Indeed, as we show 
below, the Democratic revolution of the 1930s can be interpreted in this 
fashion by means of the first-order difference equation model. 

An equilibrium value for a linear difference equation may be thought 
of asa steady state that is independent of time (May, 1974), anda difference 
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equation sequence is in equilibrium if the sequence remains constant over 
time. The discrete-time system is said to be stable if the generated sequence 
returns toward equilibrium after a disturbance moves it away from 
equilibrium. Thus, in terms of political mobilization, a d 
process is approximately at equilibrium when the support of a politica 
party remains constant through successive elections. ^ У 

The reader should understand that equilibrium does not imply stasis 
or a lack of dynamics—it does not require that all supporters stay sup- 
porters nor that all nonsupporters remain unrecruited. Indeed, so long as 
recruitment and defection rates differ from zero, defections and recruit- 
ment must occur, even at equilibrium. Rather, in equilibrium all motion 
is balanced such that defections are perfectly offset by new recruitment. 
Under equilibrium it is the ner change that is zero, not all change in the 


system. In formal terms, the number describing successive states of the 
system does not change. 


By convention, if a dynamic equation specifies some system state 
the equilibrium value for S, is denoted as S*. Note that the time subscript 
disappears, since the equilibrium value is not (by definition) time de- 


pendent. The equilibrium value for our mobilization model may be found 
in a number of ways. Perhaps the 


librium, net change (here, AM.) is z 
equation to zero, and substitute M 
solve for M* (from equation 2.7): 


most intuitive is to recall that at equim 
ero. If we then set the original dynamic 
* for М, on the right-hand side, we may 


AM. = g(L - М) - fM, 


and at equilibrium 


AM, =0 (by definition) 
so. 
0= g(L- M*)- fM* 
or 
gL g 
М* = — = —— 
P (тї). [2.24] 
And, if the first-order equation is written in standard form (from equa- 
tion 2.8): 
Misi = ao + ai M, 


it is easily shown that 


M* = ao/(1 - ai) = (gL)/(f + В) (see equations [2.25] 


2.12 and 2.13) 
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Equations 2.24 and 2.25 reveal an important feature of equilibrium. 
The equilibrium is independent of initial conditions: The process has an 
outcome that is divorced from its origins! 

As equation 2.24 indicates, the equilibrium also has a direct inter- 
pretation in terms of model parameters. Recall that L is expressed in the 
same metric as М, whereas g and f are rates or probabilities. The equi- 
librium is the product of (1) the available population, L, and (2) the ratio 
of the recruitment rate, g, to the sum of the recruitment and defection 
rates, g + f. Thus, if L=0 or if g = 0 the equilibrium will equal zero because 
no one is available to be recruited in the first instance, and because none 
of those available are recruited in the second instance. If f 7 0, the equi- 
librium will equal L because continually occurring recruitments are never 
offset by defections! Thus, equation 2.24 emphasizes the differences in 
the substantive status of model parameters: g and f are the structural 
parameters underlying the process, whereas L is a description of the 
population—a population parameter. 

The solution to the dynamic equation not only ascertains the exact 
path of a mobilization process; it also provides a basic understanding of 
stability and qualitative behavior over time. In first-order linear difference 
equations with constant coefficients, the parameter a; is crucial to this 
understanding. Depending upon the value of ai, a sequence may be 
produced that is (1) unstable and diverges toward positive and (or) negative 
infinity either monotonically or with period two oscillation, or (2) stable 
and converges toward an equilibrium point, M*, either monotonically or 
with period two oscillation, or (3) constant and maintains its initial condi- 
tion through time. These are the only qualitative behaviors that are 
possible. 

General criteria for these behaviors are obtained directly from the 
general solution for first-order linear difference equations with constant 
coefficients. Recall that this solution is given by: 


do = аа " ^ 
М = ———— +а М (fora, * 1) [2.26] 


An inspection of this equation shows that the absolute value of the 
generated series will approach infinity whenever a; is greater than 1 or less 
than -1. Conversely, if a; lies within the (-1,+1) interval, then both of the 
а? terms approach 0 in the limit. Thus, for successively larger values of n, 
the sequence approaches the limiting value, 


ao/(1 - а) = M* (see equation 2.25) [2-27] 


A constant sequence is produced by a special set of circumstances. 
When a; does not equal 1, and when 


Mo = ao/(1 - ai) = M* [2.28] 


the resulting sequence begins in equilibrium and remains so. 
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What can be said regarding the behavior of the sequence as it seen es 
toward equilibrium or diverges toward infinity? Returning to the ош 
for equation 2.22, notice that af oscillates between negative and pos E 
values if a; is negative. Alternatively, a? grows or decay аре i3 
(1.е., ever increasing or decreasing) whenever a; is positive. Thus. у 
stability (convergence or divergence with respect to рио 
qualitative behavior (monotonic or oscillatory trajectories) depen ed 
clusively on the magnitude of ai, and are not affected by either the До, | 
condition, Mo, or by ао —1ће exogenous input to the process. These results 


i icity, г scillation are 
concerning convergence, divergence, monotonicity, and oscillation 
summarized and illustrated in Figure 2.1. 


Substantive Applications 


The mathematical properties of first-order linear difference equations 
with constant coefficients have important implications for our model ol 
mobilization. We have derived the equilibrium mobilization level, M 2 
and have given it direct substantive interpretation. But is the equilibrium 
point a stable one? What are the political conditions that produce an 
unstable or divergent mobilization process? A stable one? D 

Before answering these questions, we must review the descriptive 
constraints upon model parameters. Because each parameter represents 
a substantive concept that is tied to the logic of the model, each is subject 
to certain definitional contraints. As was discussed above, both the 
recruitment and defection parameters (g and f) must lie within the zero-one 
interval. Similarly, because the upper limit (L) is a proportion, it too must 
be bounded by the zero (no one) to unity (everyone) range. 


Recall that the trajectory of M, converges toward M* if and only if the 
absolute value of a, is less 


than 1. That is, for our model to be stable 
(converge on M*), we require that: 


“S(t eye) oF “Sar С} [2.29] 


Thus, given the descriptive constraint that 


0<f,g,L <+ 


isi е. The mobilizati scribes a 
process that is inherently stable. as ; Pizauon model describe 


it must be. If our model predicted that 
i Y to positive or negative infinity over 
ume, we would clearly have to rethi negative infinity с 


(1=+=в):<0 


Z monotonic oscillatory 

[o 1 | 1 | convergent convergent 
a 

gals monotonic oscillatory 

| 1 | divergent divergent 


E 


Figure 2.1 Qualitative Behavior of the First-Order Linear Difference Equation with 
Constant Coefficients (My # M*) 


mobilization will exhibit an “erratic” oscillatory convergence toward M*, 

but if the inequality is reversed, the time-path will monotonically grow 

(if Mo < M*) or decay (if Mo > M*) toward equilibrium. Put another 

way, if (1 - f) > g, then convergence is monotonic, but if (1 - f) < в Ж. 

convergence is oscillatory§,U.E. К.Ӯ,, West Senga VES 
sr ЖШ = COs ° > ^ 


fo si = ч є 


29 

What interpretation may be given to the quantity | - f? Recall that f is 
the defection rate, or the probability that a mobilized individual ceases to 
be mobilized. The quantity 1 - f, then, gives the probability that a mo- 
bilized individual remains mobilized; that is, | - f is the retention rare for 
the mobilization process. Thus, whenever the retention rate exceeds the 
recruitment rate, the mobilization process is monotonic toward M* —the 
process moves smoothly. But if the recruitment rate exceeds the rate of 
retention, the process continually overdrives the equilibrium point, and 
convergence is oscillatory. Simply, if a party recruits individuals faster 
than it retains them, the mobilization process is an erratic, oscillatory one. 
This is an intuitively satisfying result: Given our knowledge regarding 
human behavior and the stability of partisan attachments, | - f is not likely 
to be less than g except in times of extreme political upheaval. These 
results are illustrated in Figure 2.2. 


Model Estimation 


The use of discrete-time dynamic models is not limited to the mathe- 
matical simulation of social and political phenomena; they can be used 
to analyze empirically observed processes as well. In this section model 
Seumauon is demonstrated through a consideration of Democratic mobili- 
zation in Lake County, Indiana, from 1920 through 1968. This period is an 
important one in Democratic party politics including the return to 
normalcy following World War I, the Great Depression and the New Deal, 
the Eisenhower era, Vietnam, and the political turbulence of the 1960s. 
Lake County, which includes the City of Gary, is appropriate to our 
Investigation because it has historically contained large concentrations of 
the population groups upon which Democratic ascendancy has been 
based: industrial workers, blacks, and the poor. 
oP ina oe model is written in linear form, it can, at least in 
readily du tM. by ordinary least squares (OLS). Estimates are 
but some knott oa ee on the states of the dynamic equation, 
includes 1 ыша problems Tesult whenever a regression equation 
M agged endogenous variables— prior values of dependent variables 
DT Bb he independent variables (Schrodt and Ward, 1981; 
because. stateme a шо 's an important one for dynamic modeling 
A number of p JA change explain the present in terms of the past. 
1924. aad v utions to these problems have been suggested (Hibbs, 

f we make no attempt to expand on them here. The goal of the 


res ioni y о! o 

P dt сш is only to demonstrate one model estimation st ategy: 

statist cal estimates themselves should be еу ith appr i 
spicion. valuated with propria 


The model written in standard form (from equation 2 11): 


Mui = gl * (I -f-gM 


бум tt зге г 
3-2 3. x 
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Figure 2.2 Monotonic Convergence: (g = .1, f 2.3, L = .7, Mg =.6,1-[=.7) 
Oscillatory Convergence: (g = .7, f = .7, L = .8, Mo = -1, 1 - f = .3) 
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is isomorphic to the OLS statistical model 
2.32 
Y-2mo*mX;: (+ error) [2.32] 


where пц is the slope and mo the intercept of the statistical model The 
dynamic model is estimated by OLS, with the following mappings: 


i f Lake i ing for the 
Y= M+; = the proportion of Lake County adults SQUE nds 
Democratic candidate in year t +1 fort +1 = 1924. 

1928, ..., 1968 


X = Mi = the lag of Y, i.e., the proportion of adults уоп 
Democratic in year t. for = 1920, 1924... .. 1964 

то = = 0.14 (by OLS) [2.38] 

m, 7(1-f-g)- 0.62 (by OLS) [2.34] 


We have readily obtained estimates fo 
mi, but it is more difficult to obtain par: 
from these statistical estimators. For т; 
with the statistical estimat 


r the slope and intercept, mo and 
ameter estimates for f, g, and L 
any models this is a simple task, 

ors mapping directly to the model parameters. 
Unfortunately, that is not the case here: We must solve for three unknowns 
ЮЕ from only two mapping equations. Because we have fewer equa- 
tions than unknowns (i.e., insufficient information), the mapping structure 


is underdetermined and cannot be solved without additional information. 
It is worth noting that 


even if no additional restriction can be imposed 
to determine the parameters uniquely, the equilibrium value, M*, may be 
directly obtained from t 


he empirically estimated coefficients (фо and rii) 
because (from equation 2.26): 


M* = ao/(1 - aj) 
or, in current notation, 


СЕТЕ 0.14/(1- 0.62) = 37 


аке County adult 
party. This assumption 


9 equations and two unknowns and is exactly 


determined (and solvable) as follows: 


&= то 
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f=l-mo-m [2.36] 


While this procedure does provide rough estimates of f and g, itis less than 
satisfactory: We know that L is almost certainly less than 1, not only for 
Democratic mobilization in Lake County, but for most interesting social 
processes in most environments. 

Fortunately, it is possible to utilize prior knowledge about the model 
and its parameters to derive unique estimates for all three unknowns. 
Based upon our requirement that all three model parameters must lie in 
the (0,1) interval, and upon inequalities 2.33 and 2.34, we first establish 
that: 

0< т/а < 1 [2.37] 


This inequality holds because L equals mo/g and it lies in the (0,1) interval. 
Further, because g is always positive, 


0cmocg [2.38] 
As a second step we rearrange inequality 2.34 to produce 

f=l-g-m [2.39] 
And, because f also lies in the unit interval 

OR ees el p -m <g<l-m [2.40] 


Thus, from inequalities 2.38 and 2.40 we know that an upper bound for 
gis | - mi and a lower bound is mo. 

Collecting our results, g lies in an interval bounded by mo and | - mi. 
Lacking better information it is reasonable to estimate that g lies at the 
midpoint of the (mo, 1 - mi) interval. If we assume that a normal distribu- 
tion of estimates lies within this interval, then the likelihood of choosing 
an accurate estimate for g is maximized at the midpoint. Treating our 
parameter estimates in such a manner is entirely reasonable because our 
choice of time points and elections is only a sampling from a universe of 
mobilization estimates: off-year elections, state-level elections, and so 
forth. 

Assuming that g lies at the midpoint of the (mo, 1 - mi) interval produces 
a third equation that exactly determines the system mappings. We obtain 
as an estimate of g the midpoint given by: 


g= mo + (1 - mi - mo)/2 = (1 + m - m)/2 [2.41] 


The reader should verify that our remaining mappings, when this 
estimate of g is utilized, become: 


f-I-m-gz(l-m -m)/2 [2.42] 


L= mo/g = 2mo/(1 + mo - mi) [2.43] 


Using this additional constraint, the OLS estimates of mo and mi 
produce a set of parameter estimates for Lake County that are quite 
plausible: g = 0.26, f = 0.12, and L = 0.54. Thus, the estimated Democratic 
recruitment rate during the period was twice the estimated defection rate, 
but nearly half of the eligible population was not susceptible to Demo- 
cratic mobilization.* 

Figure 2.3 displays three mobilization time series for Lake County: 
(1) the observed sequence, (2) the completely deterministic time path 
generated by the estimated parameters with the observed 1920 mobili- 
zation level as the initial condition, and (3) the time path generated by the 
estimated parameters using each observed mobilization level, Mı, to 
generate a predicted level, Mii. The estimated mobilization model 
generates a monotonically increasing Sequence that converges on an 
equilibrium value (M*) of 0.37. Notice, however, that the observed mobili- 
zation levels shown in Part A of Figure 2.3 tend to bounce around the 
predicted time path shown in Part B of the same figure 
time path that is generated using the previously estim 
the observed 1920 mobilization level as an initial condition. This is fully 
expected because we have attempted to model only the most basic mobili- 
zation mechanism without regard to short-term exogenous factors. The 
predicted time path may be thought of as an idealized path toward 
equilibrium under the condition that no short-run forces operate after 
1920. By contrast, the long-term behavior of the process including short- 


run forces, but still with Convergence toward equilibrium, is graphed in 
Part c of Figure 2.3. That picture illustrates motion toward equilibrium 
that includes the first pe 


Ч riod recoveries from successive short-term 
perturbations, Because the estimates for model parameters are based upon 
the entire set of initial conditions, i.e., each Observed point rather than 
Just 1920, Part C presents the theoretically predicted time path arising 
from the estimated model. 


a deterministic 
ated parameters with 


parameter values, On the other hand, 
à be considered an (arbitrary) initial condition. 
The model specifies the laws of behavior for mobilization in the short run. 
but the general theory of linear difference equations discloses that the 
Tun consequence, namely, the limiting or equi- 
ich the process has a tendency to move after 
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A. Observed mobilization levels in presidential elections 
Кс. 


M, 


B. Difference equation path with Mg M1920 
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Figure 2.3 Democratic Mobilization in Lake County (1920-1968) 
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ati xpressed the 
in determining the dynamics of the system. Equation 2.22 exp 
solution as 


nM a, # 
M, = ———— +а Mo (Yor a, #1) 


[2.44] 
And, on the basis of equation 2.25, we substitute to obtain 
E m Qn 
M, -M*(I =a) +а Му 
- Qn * 
LS M* +a (Mo М) 
245 
Ma = equilibrium + (decay)(initial displacement) [2.45] 


Notice that the right-hand side of 
Posing some value of the process, 
combination of the equilibrium, pl 


displacement from equilibrium measured by the quantity (Mo = M*). 
Which in turn is weighted by the п" power of the parameter ay. But ex- 


ogenous, short-run disturbances do not allow the process to approach 
gle initial condition. Rather, a whole 
» and as the effect of an “old” Mo fades 
Ws small), new initial conditions are 


equation 2.45 can be read as decom- 
say the n", into an additive ШЕП) 
из à contribution due to the initia 
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dependent; and the world almost certainly is not linear for all processes 
and all times. 

Nevertheless, the model serves to formalize a very basic social process, 
and increasingly complex models may be developed from this elementary 
building block. Increasing the complexity of a model does not guarantee 
that we will increase substantive understanding, however. Intentionally 
reduced systems models are often more insightful than their more 
"realistic" and hence more complex versions. If the modeler is sensitive 
to the assumptions they embody, first-order linear difference equations 
with constant coefficients can be powerful tools for the analysis of social 
and political change. With this caveat in mind, Chapter 3 once again 
addresses the topic of mobilization, but attempts to achieve greater 
realism by developing a nonlinear model of the mobilization process that 
incorporates the notion of social diffusion. 


3. FIRST-ORDER NONLINEAR DIFFERENCE 
EQUATIONS: PROCESSES OF SOCIAL DIFFUSION 


The social diffusion model developed in this chapter introduces a more 
sophisticated representation of change: the nonlinear first-order difference 
equation with constant coefficients. While the model is not especially 
complex either substantively or mathematically, it produces remarkably 
rich patterns of behavior over time. The model underscores the differences 
between nonlinear and linear dynamic processes. Solutions to nonlinear 
difference equations are rarely known, and no general technique exists 
for their discovery. The recursive form of this model, however, has been 
carefully analyzed, and it is possible to deduce its equilibria and qualitative 
behavior over time. (Recursive forms and other difference equation forms 
are defined in the Appendix.) n 

Models of this form have been used to explore such diverse applications 
as political mobilization, the adoption of birth control practices, and the 
spread of policy innovations. Typical applications include: (1) population 
demography, (2) social contagion and the spread of contagious diseases, 
(3) rumor spread, and (4) rapid changes in public opinion. 


Mobilization Through Social Interaction 


In Chapter 2 mobilization was modeled as a series of gains and losses 
at fixed rates. The source of change or system input, then, is constant 
through time. This constant source might be, for example, exposure to 
the media (if M, represents support of the Vietnam War) or exposure to 
air pollution (if М, represents lung cancer). But the critical assumption is 
that change is the result of constant source effects acting on individuals 
who are behaving independently of all other individuals in the population. 
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ion is enable: The 

There are many instances when this DE is птен 

underlying dynamic structure results from the RA ae 

not from independent responses to constant Source i ое 
individuals are exposed to other individuals, and this interac g 


i e dynamic 
change over time. When this is the case, the logical form of the dy 
model becomes: 


i sses * 
change 7 constant source Bains — constant source losses 
interaction gains (or losses) 


i 1 i s have been 

Dynamic mechanisms that incorporate interaction effects iod n 

demonstrated in a wide range of substantive applications. е 
suicide rates (Durkheim, 1897), voting, turnout, and partisanship( 


1927; Huckfeldt, 1979, 1980). Social interaction processes have also been 


i i ati atient 
demonstrated in such areas of behavior as choice of occupation, p 


; ; 3 ents of 
behavior in hospital wards, and racial attitudes, Dynamic ет 
Social interaction. processes have been provided by Przeworski : 
Sprague (1978), Coleman (1964), and McPhee (1963). 


A Diffusion Model 


In this section we develop a simple model of diffusion. Diffusion in 


i i RAE a constant 
human populations includes two central mechanisms: (1) a const 
Source of "exposure" 


s TUR S is 
and (2) spread from "interaction." Models of ie 
kind can be interpreted either as probabilistic/stochastic or as dete 
ministic. A model si 


t TEM t 
milar to the deterministic one developed here, put 
thoroughly Stochastic, is known in the literature as a “pure birth process 


model (Bartholomew, 1967). We are modeling change with a mobilization 
i € interpretation mi 


©. Suppose that gains are 
"social" interaction, Fina 
haracterized as 
miting, “ 
memory 


Obtained both froma constant 
lly, assume that losses are not 
"spontaneous recovery." (For 
old” rumors are no longer of interest 
erodes from simple time decay of infor- 


Source and from * 


example, disease is self-li 
and hence forgotten, or 
mation.) Thus, write 


AM, = source gains + interaction gains - decay losses p.u 


We partition the population, initially, into 
or know the information, or havi 
have not been €xposed to the di 
failed to learn the information, 


(1) those who have the disease. 
€ heard the rumor, ana (2) those who 
Sease, or have not heard the rumor, ОГ 


Thus, M, is a measure in a population of 
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those possessing the property of interest. Let Mi be measured by counting 


and express it as a proportion: 


Ms (number of individuals with the property) 
M,= — 


(total number of individuals in population) 


Similarly, those without the property may be measured as (1 - М.). Losses 
are treated as they were in the mobilization model in Chapter 2, but they 
could be elaborated in various ways. For now they are simply represented 
as -fM,, where f represents the probability of members who no longer 
have the property leaving the group. In terms of our model then we have 


AM, = source gains + interaction gains - fM, [3.2] 


Next we deal with gains. Define a quantity, L, in the metric of М, as the 
upper limit of susceptibles. For example, there are always some in the 
population who are immune and will not catch a particular disease, some 
who are not in any communication network and will not hear the rumor, 
and some who will never form an unfavorable attitude about the president. 

When Truman was at his lowest point of popularity, and when Nixon's 
popularity reached its lowest point after Watergate, the polls still showed 
that approximately 25% of the population still approved of both presi- 
dents. Thus, L is to be :hought of as an upper bound, or natural limit on 
the "spread" of М,. In this example, we think of M, as (1) starting low 
(that is, Mio approaches 0) and (2) growing in such a way that L is not 
exceeded. This second interpretative assumption furnishes a constraint 
on the parameters included in the model. Pictorially, the time path we 
have described is that displayed in Figure 3.1. Furthermore, we focus on 
“epochs” for diffusion during which conditions do not greatly change. 
Thus, we Six parameters as characteristic of a time interval of substantive 
interest, e.g., the length of an epidemic, the duration of a political cam- 
paign, or the time lapse between new and different scandals. —— 

With these interpretations, then, L defines the whole population less 


nonsusceptibles, L = | - W. And, potential possessors may be written as 
L-M 
where 


M <L always 
So we write, assuming fixed gains at rate g, 
AM. = g (L - M) + interaction gains - fM, [3.3] 


Possible interpretations of g(L - M) as gains from a constant source are 
that they could come from (1) radio broadcasts, (2) political pamphlets, 


TIME 
Figure 3.1 


Time Path for a Diffusion Process d 
R а г ote 
(3) newspapers, or (4) a contaminated well, Some readers may have n 
that equation 3.3 i 


i se 
ncludes a potential problem of double counting tho: 
who are recruited. We deal with this subsequently. 

The theoretical pr 


" izing an 
oblem, then, is that of abstractly betas ы а 
interaction mechanism in a convenient mathematical representation. 


P БАЕ А пагони air. 
est possible social interaction, the minimal ра 
Actually, most com 


munication or social interaction 
considering interactions 


NE А :mnlification of 
as only pairwise is not a serious simplification 

social reality. The descri 

a model, however, 


have the property of 
interest, and let I; denote he {а 

Similarly, let the absence 
denoted as Jo and Ji, 
interactions between | 
for example, (l, 
In matrix form: 


Property in individual J be 
respectively. Then there are four possible pairwise 
and J: (lo, Jo), (lo, J1), (ido), and (1.15), such that. 
Ji) denotes the interaction of two mobilized individuals. 


33 


Individual J 


Jo Ji 
1 + 
lo (19,49) (10,1) 
Individual I 
ij (1,9) (yd) 


This construction underlies the assumption we make about mixing in the 
Population, an assumption that is almost certainly oversimplified, but 
which is quite useful as a first approximation in representing diffusion 
phenomena. We assume, specifically, that each interaction is pairwise and 
occurs randomly and independently among individuals in the population. 
Thus, the probability that any interaction has components of a particular 
type is proportional to the frequency of the types in the population, i.e., 
Proportional to the frequency of "0" and “1” type individuals. 

In our case, one is or is not a member, one knows or does not know the 
information, one has or does not have the disease, and so forth. Clearly, 
the probability of an individual chosen at random at time t having the 
Property of interest is given by Mi. The probability of not having the 
property is | - M. But we are interested only in those interactions among 
persons who have the property and those who might acquire it through 
interaction. The expectation for those susceptible individuals in the 
Population is, again, L - My. Assuming, then, that each interaction is 
random and independent, the probability of the pairwise interaction 
(1,2) being of a particular type is given by: 


5 
EQ) = Mt 


Jo 
Kg) = L-Mt 


1 
0 
Kg) = L- Mi E(lgJg) -«L-MO(L-MO — Eoi) = м. = Mp 
L 


: - Mp MD 
FU) = My FECI, Mg) = Me М0) ка.) = MoM 
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The construction given above supplies our basic interactive rule. 


Interactions of (M, M;): no spread, 

Interactions of [(L - M),(L - M,)]: no spread, 

Interactions of [M,,(L - My]: yield spread possibility, and 

Interactions of [(L - M), Mi]: yield spread possibility. 
From the matrix above it is evident that the relative frequency of pairwise 
interactions in the Population between a mobilized individual and a 
potential recruit is given by: 2(L- M)M.. It is more convenient, however, 
to drop the doubling of the (L - M;)M, quantity, and force the coi rection 
Onto the parameter that gives the rate at which these "across group 


interactions cause movement into or out of M, We denote this rate as the 
parameter s, and write: 


AM: = &(L - M) + sM(L - М) - fM, [3.4] 


It would appear that we are finished. N 
equation 3.4 might result in an individ 
constant source effects and from interaction effects, Since an individual 
can become an M, only once from time t to t + l, we must correct for the 


potential double counting of gains from both constant source and inter- 
action mechanisms, 


otice, however, that the form of 
ual being mobilized both from 


action, (L - Mj), remove those already 


recruited from constant source 
effects in an instant of time. Hence, writ 


e 


[KL - M)- g(L - Mj] 


as the pool available for Social interaction. 
Our final elaboration of the diffusion model may be written, then, as: 


АМ, = g(L - M) + sM((L - M) - g(L - М))] - fM, [3.6] 
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Slightly more complex versions of the model are possible if we elaborate 
the loss term -fM,, but we shall continue our analysis by considering 
only the most elementary case (from equation 3.6): 


AM, = g(L - M) + sM[(L - M) - g(L - MJ] - fM. 
Algebraic manipulation yields a quadratic form for the model: 

AM, = (sg - 5) М? + (sL - f- g - sgL)M: + gL [3.7] 
or by disaggregating the AM, into Ми: - Mi, obtain 

Mia = (sg - 5) М2 + (1 + sL - f- g- sgD)M + gL [3.8] 


The reader should verify that equation 3.8 is a first-order difference 
equation with constant coefficients. Note, however, that the dynamic 
equation is nonlinear, since one of the system states is squared (i.e., 
second-degree). The model, then, is a first-order quadratic difference 
equation with constant coefficients. 


Parameter Restrictions 


Descriptive adequacy imposes at least the following restriction: 


[3.9] 


may be required. From the model, 


but stronger restrictions for g and s < ! 
the following inequality: 


the total gains evidently must satisfy 


0< g(L - M) + sMi[(L - М) - g(L - Mo] (1.5 М) [3-10] 


for: жь апа 


0« 
0cg«l 


And by imposing positive s: 0 < $ [3.11] 
1 gains (change) must be positive 


Inequality 3.10 simply states that the tota be positi 
4 S pool of susceptibles. Simplifying 


or zero and less than or equal to the entire 
Inequality 3.10 we have 


0<ре+5М(1- 8) <! М. #1. [3.12] 
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which may be written 


0< g(1 - sM) € (1 - sM.) yielding [3.13] 


0<р<1 [3.14] 


Which is exactly the same restriction assumed earlier. Similarly, we can 
solve inequality 3.10 for the restriction on s. Thus, we have 


OS g/(1- g) + s < 1/(1 - g) for g not equal to 1 [3.15] 
Which reduces to 
0<5<1/М, with M.>0ands>0 [3.16] 


If s < 0, then inequalities 3.11 and 


3.16 need no longer hold because the 
result of social interaction is to г 


emove members from the population who 
have the property of interest, and the interaction term is a loss term. The 
reader might wish to develop the analogous parameter restrictions for the 
case s < 0. 


Model Estimation 


This model turns out to be nonlinear in its parameters (f, g, L, and s), 
but estimates for the parameters can be obtained if we impose 
hypothesis about L. The easiest solution to the Problem is to 
to I. That is, make the entire Population susceptible for recruit 
set L at a particular value, the model is stil 
now, unique estimates can be found. 


algebraic manipulation, we can put th 
following form: 


an a priori 
set L equal 
ment. If we 
1 nonlinear in parameters but 
If we disaggregate AM, and use 
е model in equation 3.7 into the 


Mei = Gg - Mi + (1+ sL- f- g- sgL)M, + gL [3.17] 
This is equivalent to the least Squares estimating form 


Y = m2X2 + mXi + mo (+ error) 


[3.18] 
where: 
Y =Ma 
X:= Mi 
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It is instructive to note that the only set of observations required for 
estimation of this model is a time series of one state variable plus its lag. 
Using a time series of a single political variable, for example, we can 
investigate the relative contributions from constant sources (g and f) and 
from social interaction effects (s) with the same model. Least squares 
procedures may be used to obtain estimates of the coefficients m; in equa- 
tion 3.18. The estimating relationships for the parameters are specified 
by the system 


sg-s -m [3.19] 
I*sL-f-g-sgL-m 
gL = mo 
For the simplest case assume that the entire population is susceptible, then 
L= I. With this assumption system, equation 3.19 can be solved uniquely 
for the parameters s, g, and f in terms of the coefficient т; as set forth below 
f=1-m-m-m [3.20] 
g= mo 
5 = m/(m - 1) 


Even though it is relatively easy to obtain values for the parameters 


using least squares estimation procedures, an explicit solution for the 
quadratic form of the first-order difference equation is not available. 
Conditions for convergence or divergence are known, however, às are tests 
that determine the qualitative nature of the time path. Furthermore, if 
the process is convergent, an explicit expression for the limit is available. 


An Empirical Example 


For illustrative purposes we have estimated the model using a time 


Series of voting data for Essex County. Massachusetts. The dong 
exami i the vote used is Democratic vote 10 s 
edere pes of the presidential 


dent, with M, defined as the Democratic proportion c 
vote. (Data were collected from America At The Polls and America Votes, 


edited by Richard Scammon.) This covers the same period of aon 
mobilization as that used in Chapter 2, but uses a northeastern Democra 
suburb of Boston. 
Essex County parame 
efficients using equation 3.20, are presente 
Least squares coefficients: 


ter estimates, calculated from least squares CO- 


d below. 


4 [321] 
tho = 0.16 5-095 th = -0.45 
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Parameter estimates: 


whenL-l, then f=0.34 g=0.16 and 5= 0.54 
when L=.74, then f=0.25 g=0.22 and s-0.57 [3.22] 


Two sets of estimates were calculated from the m, 


and .74, respectively. The value .74 for L might be 
residual support for Truman and N 
Furthermore, 74% 
Essex County durin 


by setting L equal to | 
justified by noting the 
ixon at their lowest was about 25%. 
is the highest Democratic Proportion of the vote in 
g this period. Thus, it is perhaps reasonable to suppose 
about one-fourth of the Voting population will never move to any cause— 
in Essex County, they will never become Democrats. The parameter 
values produced when L is set equal to one indicate that social interaction 
effects are Stronger than individual level effects. (Essex County has a 
strong Democratic organization, which suggests that gains by the Demo- 
cratic party are a result of interactions between Democratic supporters 
and recruitable nonsupporters.) The individual-level effects are also 
important, Suggesting that national-level forces—social, economic, and 
vidual voters in the county. Thus, 
ual level change at the local level. 
: mates when L is set to .74. Social 
interaction effects аге į ightly, but individual loss effects are 
lier studies (Berelson et al., 1954; 
ated the importance of personal 


Priori restrictions ( 
M 
analyzing the mobilizatio 
ам re es n aluated using appropriate 

ution due to the presence of lagge: iables 
(see Hibbs, 1974), "PAROI varla 


as they do in this instance), 


Analyzing the Model 
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The problem is compounded by the fact that, unlike the linear dynamic 
system considered in the first chapter, nonlinear systems frequently ex- 
hibit multiple equilibria rather than a single, globally stable equilibrium 
point. Furthermore, the best general strategy is to ascertain the dynamic 
stability of some "neighborhood," or local region, around the equilibrium 
points. Only if the dynamic equations are linear will neighborhood stability 
and global stability necessarily be identical. This means that, whereas 
initial conditions do not affect long-run qualitative behavior in linear 
dynamic systems, they can critically alter long-run behavior in nonlinear 
dynamic systems. 

Another complication that may arise from nonlinearity in the dynamic 
equations is that equilibrium need not be a point, but can be a stable 
limit cycle. As the term suggests, in a limit cycle the states of the system 
undergo cyclical patterns of change across time, and a stable limit cycle 
will return to its oscillatory trajectory after local perturbations (for more 
detailed technical discussion of these nonlinear system characteristics, 
see Rosen, 1970; Hirsch and Smale, 1974; and May, 1974). 


Equilibrium Analysis 


its states is zero over 


If a process is at point equilibrium, net change in 
will be the same. 


time, Thus, at equilibrium, all successive values of Mi 
That is, for a dynamic equation of order k, at equilibrium 


Mux = М, = M* [3.23] 


tationary point). Wecan 
23 by setting AM: equal 
nd then solving 


where M* is an equilibrium point (also termed a si 
then determine all values that satisfy equation 3. 
to 0, substituting M* for all M, in the dynamic equation, a 
for M*. Utilizing this procedure for our model gives 


0 = (sp - s)M*? + (sL - f - g- sgDM* + gL [3.24] 


Equation 3/24 is simply a quadratic equation in M*, which may be 
solved by the quadratic formula. Recall that for an equation 


0 = ах? + ЪХ tc [3.25] 
the solution for X is given by 
X = (-b + Vb - 4ac )/2a [3.26] 
Equation 3.24 maps directly to 3.26, with 
[327] 


а=5р- 5 


40 à; 
3.28 
b-sL-f-g-sgL [ 
L [3.29] 
cig 


i а lex 
Thus solutions for equation 3.24 are given by the rather comp 
equation: 


M* -[-GL - f- g- set) 


+ VGL-f-g-sgL) -asg aL 1/2(sg - s) [3.30] 


М» = [(.05) + №005) - 4(-.45) (16) ]/2(-.45) [3.31] 


Which reduces to 


М* = [05 + 2005 ]/C.9) 
M* = 0,54, -0.66 


Same, i.e., the proces: 


S remains dynamic but the 
à number— no longer 


changes.) 


Local Stability Analysis 


A logical next question is; What happens when the process is disturbed 
by external forces? These may include Short-run Political forces such as 
political scandals, economic conditions, and international crises. These 
short-run forces would shock t 


he system away from the equilibrium. After 
these forces have perturbed the political System, will the process converge 


toward the equilibrium Point or will it diverge? We have solved the 


quadratic form to determine its equilibria, but we have not yet provided 
a means to test for stability, 
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For nonlinear dynamic equations stability can be local or global. 
Local stability means that, within some neighborhood of the equilibrium 
point, the process will converge toward the equilibrium point after it is 
perturbed. Global stability implies that the process is stable no matter 
what the perturbation (for a more complete description of stability see 
Rosen, 1970; May, 1974). 

First we apply a technique to investigate local stability at theequilibrium 
points, and then we will take advantage of some general results of Chaundy 
and Phillips (1936), reworked by Sprague (1969), to determine global 
stability for this specific quadratic form. The technique used to investigate 
local stability has a wider application for more complex, nonlinear 
models. 

^n analysis of small perturbations around the equilibrium point M* 
begins by writing the perturbed mobilization level as 

AM, = M* +X, [3.33] 
Here X, measures a small disturbance to the equilibrium M* within some 
Specified neighborhood. Now X, must obey the original laws of change for 
the mobilization process given by equation 3.6. However, in a small 
neighborhood of the equilibrium—X, is a small perturbation—a basic 
result from the calculus of numerical analysis may be invoked. It is possible 
to represent the process by a particular transformation, and, with an 
additional trick justified by the assumption that the perturbation 15 small, 
characterize the behavior of the nonlinear equation as essentially linear 
in the neighborhood of the particular equilibrium point chosen. This 
Procedure is called an expansion in a Taylor series. | . 

An approximate difference equation for the perturbation measure is 
Obtained by a Taylor expansion of the equation for our original model, 
equation 3.6, about the equilibrium point. The Taylor series expansion 
Provides a linearization of the model because the linear term in the 
expansion dominates the series in a small neighborhood, and terms of 
ш two and higher сап be neglected. The expansion takes the following 
orm 


3.34 

АХ‹=аХ, or Xac(*taX pou 

Where a is the partial derivative of AM, with respect to M, evaluated at the 
equilibrium point M* (obtained from equation 3.24). Hence. 

[3.35] 


a = [Q(AM)]/(0M) = 2(sg - 9M* + sL - f- g- 581. 


; i i ighborhood 
It measures the mobilization growth rate 1n the immediate neighb 


of the equilibrium point. i i 
Equation 3.34 is a first-order linear difference equation for which we 
have an explicit solution. It has the form 
[3.36] 


Xı = ХАІ + a)! 
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where Xo is the initial small perturbation. The disturbance dies away if 
(1 + a) lies in the open interval (-1,+1), diverges if (1 +a)> 1ог(1 +а)<-1. 
and is constant if (1 + a) = 1. Thus the neighborhood stability analysis of 
the equilibrium point M* shows the point to be stable if and only -1 < 
(1 + a) € +1, or more simply, -2 < a < 0. 

Applying these results to our empirical example, the mobilization of 
the Democratic party in Essex County, substitute equation 3.35 into 
equation 3.34 to obtain 


AX, = [2(sg - s)M* + sL- f - g- sgL]X. [3.37] 
Disaggregating AX, gives 


Хи = [1 + 2Gg - s)M* + sL - f - g - sgL]X. [3.38] 
where 
1 +a=[I + 2(sg-s)M* *sL-f- g- spl] [3.39] 


Evaluating the coefficient of X. at M* = 0.54 using the estimates for the 
parameters s, f, g, and L from equation 3.22 gives 


(1 + a) = (1 ~ .49 - .05) = 0.46 [3.40] 


Since l+a 5 -46 and therefore lies in the (0.1) interval, we know that the 
disturbance is monotonically convergent and dies out. This means that 


the equilibrium M* = 0.54 for Essex County i 
і nt ili 
zation of the Democratic party conver; Re i nav инан: 


43 


Bence, and ultimate qualitative behavior can be developed from their 
discussion. Only a few of the results are presented here; the inquisitive 
reader should consult the original source. 

First define a quantity K by 


K =[-1 + V1 + 4((B/2) - B/2 - AC]]/-2 [3.42] 


Where A, B, and C are from equation 3.41. This produces three possibilities: 
two real and unequal Ks, two real and equal Ks, or a pair of complex Ks. 
Six cases are considered below. 


Case 1. If K given by equation 3.42 is complex, then the process is 
divergent, diverging to infinity. hi 

If K is real, choose that K that is > .5. One of the Ks should meet this 
condition. 


Case 2. If | AMo + B;2| > К, then the process М, diverges to infinity. 


Case 3. If | AMo + B/2| = K, then the process M: is stationary. This does 
not mean the process will converge if displaced. 


Case 4. If |AMo + B/2| < К and 1/2 € К € 3/2, then the process Mi 
converges to a value 
M* = (1 - K - B/2/A [3.43] 
The limit M* is dependent on A, B, and C since K depends on C. Con- 
vergence in this case is monotonic if 1/2 € K < 1. 


Case 5. If |AMo + B/2| < K and 3/2 < K <2, then the process М, oscil- 
lates finitely. 


Case б. If |AMo + B/2| < K and K > 2, then the process M: diverges to 
infinity with a certain exception, i.e., if Mo is chosen so that the роге 
[АМ + B/2| is ап element of a set involving the square roots eae tae 
sion К“ - K, then the process M. oscillates finitely (this discussion E base 
on results from Sprague, 1969). Given particular estimates ofany quadratic 
form, we may deduce its qualitative behavior over time. 

We now apply these conditions for convergence to the data on Demo- 
cratic mobilization for Essex County. We have already determined that 
there is a locally stable equilibrium at M* = .54. We now take advantage 
of the preceding results to see if the locally stable equilibrium satisfies the 
Chaundy/ Phillips conditions for global stability. First the real numbers 
A, B, and C are defined 


А = п = sg-s 
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B-m-l*sL-f-g-sgL 
C=m=gL [3.44] 


Substituting the estimates for the parameters for s, f, g, and L for Essex 


County from equation 3.22 into the formulas in equation 3.44, we obtain 
the values А =-.45, B = .95, and C = 


-16. Using these values we calculate K 
as follows 


K =[-1 x V1«4(23 - 475 + 07) ]/-2 


K = 78, 23 [3.45] 
We choose .78 as the value for K and find 


amine the value | AMo + В/2]. Substituting 
|.45Mo 


that Case 4 applies. Now ex- 
values for A and B we obtain 
+ .475|. The condition for convergence of the process is 


|-.45Mo + .475| < 78 [3.46] 
The inequality 3.46 emphasizes that convergence 
equations is dependent upon initial conditions. T 
the mobilization process is an important consi 
nation of long run limitin 
equality 3.46 hold? We 
Here, if inequalit 
all values of Mo. 
values 0 and 1 for 


Jor nonlinear difference 
hus the starting point of 
deration in the determi- 


9r convergence. We also note that con- 
Vergence of the process i ic si ies i ,1) interval. 
Finally, M* calculated using equation 3.43 equals .54 for Essex County. 
This is the same value obtained using the quadratic formula, which is as 
it should be, 


Summary 


We have investigated the m 
quadrati 
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systems. This, of course, is not the case with linear models, since neighbor- 
hood analyses describe global dynamics as well. In nonlinear systems, 
however, we do not know what will happen if the system is disturbed with 
sufficient force: The states may return toward equilibrium, move into a 
stable limit cycle, or explode to infinity. The qualitative behavior will, 
in fact, depend upon the second- and higher-order terms of the Taylor 
expansion, which are discarded in a neighborhood analysis. Fortunately, 
dra case of the quadratic, some general results are known. Conditions 
onvergence and divergence were reported and are applicable for all 
models that take this form. In this particular instance, then, a partial 
global analysis is possible. 
We continue our treatment of dynamic models in Chapter 4 by returning 
to linear difference equations, but we extend our consideration to systems 


of interdependent equations. 


R SYSTEMS AND HIGHER-ORDER EQUATIONS: 
ARMS RACES AND FEEDBACK PROCESSES 


In previous chapters our development of discrete-time models is limited 
19 single first-order equations. This chapter extends our consideration to 
higher-order models. Two different models are developed: the Richardson 
arms race model and a feedback model of legislative review. The two- 
nation, linear arms race model can be reduced, for analytic purposes, to 
a linear, second-order difference equation with constant coefficients, and 
the feedback model treats single-equation, higher-order lag processes 
More generally. During the course of this chapter a matrix representation 
of the arms race model is presented, some useful analytic approaches are 
explored, and general solution strategies for linear models or linear 
systems of higher order are outlined. 


Higher-Order Models 

er it is dynamic or static, mathematical 
or verbal, serves as an abstraction of reality that necessarily simplifies the 
Subject of study. A first-order dynamic model simplifies time dependence 
by expressing a current system state as à simple or complex function of 
the immediately preceding system state, ignoring that state’s extended 
history. The reality is clearly more complex than the abstraction. Given a 
System state at t, we might reasonably extend the functional form in- 
definitely back in time so that: 


Any social science model, wheth 


S = f(X Sex) k= 1,2, 8s [4.1] 


ped statistical technologies for identifying 


Econometricians have develo 
ese statistical techniques 


the structure of such time-ordered dependence. Th 
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are useful as an inductive exercise in model building—as a way to develop 
substantive ideas regarding the nature of time dependence and system 
memory. A theoretical model, however, as it is understood in this mono- 
graph, is properly seen asa direct mathematical expression ofa substantive 
idea. Thus, it is an inappropriate modeling strategy to extend the structure 
of temporal dependence simply to be safe, or simply because the data are 
available. The structure of time dependence ultimately must be specified 
on the basis of a theoretical understandi 
In many instances this means that a first 
favor of a higher- 


ng of the phenomenon in question. 
-order model must be discarded in 
order model. For example, in the mobilization models 
of the previous chapters it might be appropriate to argue that recently 
recruited individuals (those recruited in the last time period, t - 1) have 
different defection rates than individuals with longer membership his- 
tories (those recruited at t - 2 and before). This would imply two loss rate 
hypotheses and a second-order process. 

Higher-order models can also be developed on the basis of two or more 
separate, but functionally interdependent, statements of change. These 
component parts may be first-order specifications: Two interdependent 
first-order models produce a second-order System, three interdependent 
first-order models imply a third-order system, and so forth (Samuelson, 
1974). The Richardson arms race model developed in this chapter is written 
as two first-order equations that produce higher-order reduced forms 
(a single equation of order two for each system state). Thus, in introducing 


Second-order models the concept of linear interdependence of system 
States is also introduced, 


The Richardson Arms Race Model 


Lewis F. Richardson's mod 
Spectrum of dynamic modeling 
ments buildups (see Gillespie an 
his model as a continuous-ti 


ДХ, = kY.- aX + g 
where: 


X = the armament level in nation X at year t 
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Y; = the armament level in nation Y at year t 


a =the economic burden experienced by nation X in attempting to 
maintain its armament level 

k = the threat felt by nation X due to the armament level of nation Y 

g = the effect of ambitions, prejudices, and hostilities upon the arma- 


ment level of nation X 


In words, nation X will modify its arms level from time t to t + 1 as 
follows: It will examine the arms level of its enemy, and (disregarding 
other considerations) respond by building new arms at rate k. Nation X 
will also consider the costs of maintaining its current armaments (X), and 
reduce its arms levels at rate a on the basis of its own perceived economic 
constraints. Finally, purely domestic considerations—aggressiveness, na- 
tionalism, pacifism, patriotism—create pressures to increase (or decrease) 
armaments by an amount g independently of what nation Y is doing. 
It is the net effect of these three considerations that produces the change 
(positive or negative) in X:. | 

Simultaneously, nation Y makes similar decisions regarding its own 
future armament levels according to the same rules. With only minor 
modification, the armament behavior of nation Y can be expressed as: 


AY, = k’X, -aY +g [4.3] 


where: 

- the economic burden experienced by nation Y in attempting to 
maintain its armament level 

7 the threat felt by nation Y due 


the effect of ambitions, prejudic 
ment level of nation Y 


= 
П] 


to the armament level of nation X 
es, and hostilities upon the arma- 


Notice that k, k', a, and a' are all rates or proportions, while g and g' are 


both /evels measured in the same metric as X, and Y.. Р mem 
Clearly, nations X and Y are locked inan interdependent system delin 


by equations 4.2 and 4.3. Furthermore, this system is simultaneous in > 
timing. Interdependence and simultaneity are typical of many m 

Processes: The behavior of a first actor depends on the behavior of a 
Second actor and vice versa. Game theory is one well-known шсш 
analyzing such conflict situations, but dynamic models are also powerlu 


tools for examining them. 


Arriving at a Second-Order Reduced Form 


The two-nation model is composed of two equations (4.2 and 4.3) that 
are both first order, but the equations are clearly interdependent. eum 
Of the first-order equations can be translated into the recursive form so oe 
in Chapter 2 because both have exogenous system states on the right-han 


48 


side. In order to explore the dynamics of the system represented by these 
equations, we translate them into a different form. First, by disaggregating 


the change operator, A, and advancing equation 4.2 one time unit, we 
produce: 


Xoz = (1 - a)Xu + KY +g [4.4] 
Algebraic substitution of equation 4.3 into equation 4.4 yields: 
Xo = (17 a)Xii + k[(1 - a) Y, + k’X, + g]*eg [4.5] 


Our goal is to express Х‹ as a function of previous X. alone, but equa- 


tion 4.5 still contains Y, on the right-hand side. Rearranging equation 
4.2 produces: 


Yi = [Xii - (1 7 a)Xı - g]/k [4.6] 


Finally, by substituting equation 4.6 into e 


quation 4.5 and algebraically 
rearranging the result, we arrive at our go 


al. 
Xuz = Q7 a-a)Xi + (Kk'- 18a & a'- aa) X, + kg’ + ga’ [4.7] 

The armament behavior of nation X is now expressed solely as a func- 
tion of its own previous armament levels, even though the original theo- 
retical model clearly recognizes the interdependence between nations X 
and Y. The price of obtaining an analytic form for the behavior of X, that 
is not explicitly dependent on Y, is an increase in the order of the lag struc- 
ture for the new recursive form. That increased time lag implicitly captures 
the interdependence of the original system. With appropriate rearrange- 


ment, equations 4.2 and 4.3 can also be used to express the armament 
behavior of nation Y in a recursive form: 


Ү = (2 - a -ауү,,, + (Kk-1*a*a'- a’a)Y,+k’g-g’a [4.8] 


Putting the Arms Race in Matrix Notation 


The arms race model is easily expanded to accommodate interrelation- 
ships between more than two nations. The notation for the multination 
model quickly becomes burdensome, however. Just as a two-nation arms 
race produces a second-order System, an arms race with n nations produces 
à system of order n. A more convenient expression for such a system 
makes use of matrix algebra. 

A matrix generalization is most easily accomplished (and made 
extraordinarily useful) if it is recognized that the operation of taking a 
first difference (АХ, = Xw - X,) defines a linear operation. Thus, delta (A) 
is a linear operator, and for many purposes it may be treated algebraically 
as an ordinary constant. The theory of linear Operators is part of opera- 
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at (Churchill, 1972) and is developed for the discrete- 
linear ps depen, ( 1958) among others. Most techniques for solving 
Mis. cen en © "E difference equations rely on operational methods 
ec aba pectively, the Laplace and the so-called z-transform. In prob- 
i e же these methods arise when generating moments of a probabil- 

à ribution. All that is used here, we reiterate, is the basic fact that the 
Operator delta may be treated as an algebraic quantity. 

Rewrite equations 4.2 and 4.3 as a simultaneous system 


AX = -аХ +kY tg 


AY, =-a'Y, + k'X, +р [4.9] 


Next move the state variables (X. and Ү,) to the left-hand side, obtaining 


(А +а)х, -kY,-g 
-kK' X, + (A +а)ү sg [4.10] 


eral for the n by n case, 
sa matrix and to write 
tors. Thus, 


" Our strategy at this point, which is perfectly gen 
is to detach the coefficients of the state variables а 
the state variables and the right-hand side quantities as vec 
write 


(А +а) -k X, } g 
-k' (Ata) Y, g [4.11] 


Which expresses our system in matrix form. Pursuing our plan of treating 
the “A” as a constant we may invoke rules of matrix algebra and solve for 
System states X, and Yı. Readers familiar with matrix algebra may wish to 
Solve equation 4.11 for system states at this point. Readers less familiar 
with matrix techniques are directed to the Appendix where we find the 


Solution using Cramer's rule. 


Finding the System Equilibrium 


А Even without matrix operations, the equilibrium point in a two dimen- 
sional (X,Y) plane is readily obtained using the system in equation 4.9 
and simple algebraic manipulations. Replace AX: and AY, with 0 to denote 
à steady state, and replace X; and Y: with X* and Y* to denote the unique 


equilibrium point that produces the steady state. 
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0--aX* + kY* +g 
0 = -a'Y* + k'X* + g [4.12] 


In order to solve for X*, isolate Y* on the left-hand side of both equations 
in the system. 


-Y* = -(a/)X* + (8/0) 
Y* = (K'/a?)X* + (2а) pend 


Adding the two equations together and algebraically rearranging the result 
gives the equilibrium for X: 


х*= St tek [414] 
aa’ - kk 
Parallel procedures produce the equilibrium value for Y: 
gk + gia [4.15] 
eS : 
aa’ - kk’ 


Note that while equations 4.14 and 


4.15 give equilibrium values for X 
and Y, the system is in e 


quilibrium only when both X, and Ү, are simul- 
taneously at equilibrium. We t 


hus denote the equilibrium for the system 
as the point (X*.Y*). 


The system e 
sional 


: i will be the qualitative nature 
of the path approaching (or diverging from) equilibrium? For the case of 
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linear di x А 

sus eye equation: even without constant coefficients, solutions 
specified. The d an appropriate number of initial conditions are 
н THEE у aining issue is finding such а solution, and that is 

ask. 
I ši "Т 
tcc of higher-order solutions we rely on a general 
thar the tee ni o linear difference equations. This theorem asserts 
чай исо: sa p fora difference equation is the sum of any par- 
(Goldberg, es the general solution for the homogeneous form 
SUBE аи tte к а to apply this theorem, a recursive form differ- 
on Thé орана e iram order is rewritten with all states and state lags 
oPGElSm i x ^S e and all other terms on the right-hand side. The 
ihis мен п roken down into two parts: (1) finding any solution for 
ЖЕШСП ie ed and (2) finding the general solution for the homo- 
Didius a Fes for the complete equation with zero right-hand side. 
thé "ipaa к еза ideas we first demonstrate this solution strategy for 
viously: salve linear equation with constant coefficients that was pre- 
3 in Chapter 2. 
Consider the first-order equation in its recursive form 


Visage ai [4.16] 
This may be placed in canonical form as the complete equation 
Yu 7 arY, = ao [4.17] 
with the associated homogeneous equation 
[4.18] 


Yur -aY = 0 


know that one particular 


Fro n 
m earlier developments (see Chapter 2) we 
uilibrium solution: 


solution to the complete equation 4.17 is the eq 

Y* = ao/(1 - а) [4.19] 
quation 4.18. In order to apply the 
be obtained for this homogeneous 
the reader to Goldberg (1958) 
d that they are unique once 


N Р 
ihes. consider the homogeneous € 
fone ae a general solution must also 
faci We exhibit such a solution and refer 
d pea that such solutions always exist, an 
p conditions are supplied. 
fun € attempt to find a solution for e 
ctions and suppose that 


quation 4.18 among the simplest 


[4.20] 


may work. Substituting equation 4.20 into equation 4.18 yields 


[421] 


+ 


1 
г" -ar=0 
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which, when divided by r' (for r' # 0) gives, on rearrangement, г = ai. А 

In order to allow for different initial conditions, a constant that varies 
with initial conditions somehow must be introduced. Try a simple multi 
plicative c. If aj is a solution to equation 4.18, is it also true that cai isa 
solution? The reader should substitute and see that cai is indeed a solution. 
In fact, it is the general unique solution with c an arbitrary constant to be 
determined by initial conditions. 

According to the theorem, the solutions to equations 4.17 and 4.18 may 


be combined additively to arrive at the general complete solution. We 
try obtaining 


Yı = Y* + cai 


[4.22] 
Suppose an initial condition is Yo. Then we have 
Yo = Y* + са? [4.23] 
which implies 
c= (Yo - ү») [4.24] 
Thus equation 4.22, our solution, may be written as 
Yı = Y* + (Yo - Y*)at [4.25] 


just exhibited will work i 
frequently onerous, 


1 » our only task is t 
equation. That equation is, for state X, 


Men (а t2 2X (раа ke, a)X 2 0 [426] 
Let 


b-a*a'-2, 


С Haat ие аде 
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and 
Z = either X, or Y, 


о " s 
ur problem may then be recast as finding the general solution for 


Zu2 + bZ + cZ, = 0 [427] 


P ed as чу -— А 
Toceed as before and try г". Substitution and division by r' yields 


r+br+c=0 [4.28] 
aun 4.28 is recognized as a polynomial in r of order two—the qua- 
dratic, It is frequently referred to as the characteristic equation of the 
System; Evidently, either root of the quadratic will satisfy equation 4.27 
en before, we multiply each by a constant to allow for the two sequen- 
2 initial conditions required by the general theory. Thus, assuming the 
0015 of equation 4.28 are real and unequal, the solution to equation 4.27 
may be written as 


Zi = Z* + cur) + cord [4.29] 


Equation 4.29 may be used to determine the values of the constants given 
initial conditions Zo and Zi. The appropriate relations are 


€,7Z,-Z*-c, 


e; =(Z, - Z9)f(r, ту) - Ir (Zo - 29M 7) 


-z*)- Е 
а Ма = [4:30] 
t3 =i 


Note that the form of ua іол in equation 4.29 is similar to the form 
of the solution shown in equation 4.22 except that an extra root and 
Constant occur. This is general: One new root will be added along with the 
new constant for each increment in the order of the system. a) 
Before turning to the analytic potential of this solution, two additional 
complications must be considered. First, it may be that the solution will 
have only one (repeated) root. Second, it may be that the roots are imagi- 
nary, in which case they occur as complex conjugate pairs. For present 
Purposes note that these cases are technically manageable, leading to 
Unique solutions that exist in the real domain. If the roots are complex 
then the solution contains a trigonometric function that supplies an 
Oscillatory motion to the states of the system with period longer than two. 
_ The most important fact about these potential difficulties, however, 
is that they exhaust the possibilities. Fundamental mathematical theorems 
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provide assurance that the roots of a characteristic equation with co- 
efficients that lie in the real-number domain will produce roots that are 
either: (1) unique and real, (2) repeated and real, or(3) complex conjugates 
with a trigonometric real representation. 

Extracting the roots of characteristic equations with order greater than 
two may beatime-consuming task in numerical analysis (Hamming, 1971). 
Fortunately, the solution form can be used to obtain a set of criteria for 
determining stability expressed solely in terms of inequalities involving 
coefficients of the homogeneous equation (May, 1974; Samuelson, 1974). 
Less fortunately, if the qualitative nature of the approach toward (or 
departure from) equilibrium is to be determined, it is necessary to know 
whether the roots of the solution are real, complex, and (or) repeated. 
Knowing the motion in the system requires solving the characteristic 
equation for its roots, and this may be hard work if the system is at all 


large, requiring numerical approximation techniques and a high-speed 
computer. 


Putting the Solution to Work 


Inspection of equation 4.29 casts light on both the quality of time paths 
and convergence. Clearly, for the case of real roots, all roots must be less 


than one in absolute value if the process is to approach the limiting or 
equilibrium value Z*. If one 


the time path of the System s 
oscillate with period two. If 
single quantity can be devel 
real and imaginary parts of 
modulus of the roots in the c 


; » that the time path of th 
with period greater than two.* 


b? - 4c 


b? 
net d S [4.31] 


But, on the basis of the equalities shown under equation 4.26, we can write, 
after some simplification, 


3 


- 


Th i 
us, trigonometric motion will result if 
(а-а) + 4kk' « 0 [4.33] 


Which simplifies to 


4kk’ € -f(a - а) 
ied paca [4.34] 


In " 
ния oscillation of armament levels depends on the relation- 
опот E parameters attached to system states—a nation's own 
fidt depend ax | ion and its fear of an opponent's armament—but does 
тене далар internal grievances (the g and g terms). But the deductive 
à. and a' are ge from the form of inequality 4.34 is stronger. If k, K’, 
ane its mte à ушш than zero (which must be true if the model is to 
4:34 Scale : e interpretation), then the left-hand side of inequality 
long-period uir un while the right-hand side is ar most zero. Thus, 
model. The inre ае > impossible given the logical structure of the 
term (period As next quen s Does the model also rule out short- 
een ^ ) oscillation? Period two oscillations will arise provided 
—— 
(a+ a’- 2) - V/(a- а) + 4kk <0 [4.35] 
riod two oscillation on a priori or 
5 can be satisfied if internal eco- 
nd perceived dangers sufficiently 


нн does not preclude ре 
ТОНДО пош because inequality 4.3 
ез rains are sufficiently small а 
ub the formalism leads to de 
Gffsran i interpretations. And thes 
Йй eale "cuni means for empirica 
Bütór ini стак: measures of statistical fit 
rates be compare the qualities of real- n 
Mirada are predicted on the basis of the hypothesized m 

quate model should reproduce the qualities of an ol 
process. 


ductive conclusions that have clear 
e deductively obtained conclusions 
1 investigation and analysis. Rather 
or correspondence, the investi- 


world arms races to the arms 
odel. In short, 


bserved social 


Stability Conditions 


Goldberg (1958: 171-172) develops necessary and sufficient stability 


conditi 

ditions for the second-order difference equation: 
Condition I: 1 *b*c71 

coefficients yields 


Substituti 
bstituting the Richardson model 


sategal ead ke - gc > 
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which simplifies to 
aa’ > kk’ [4.36] 
Condition II: 1-b+c>0 

Substitution yields 
l-a-a'*2*1*aa'-kk-a-a'»0 


which simplifies to 


4 - 2a - 2a' + aa^ > kk’ [437] 
Condition Ш: 1- съ 0 


Again substitution and simplification yield 


kk’ * a * a' > аа’ [4.38] 
Now, it can be readil 


y shown that as long as the four model parameters 
(a, a’, k, k^) lie within t 


he (0,1) interval, the conditions in inequalities 4.37 
€ satisfied. We have already argued that the parameters 
‚ and two examples show that it is entirely reasonable 


i tum. In substantive terms, a balance of power 
results. If the interaction between the economic strain of the two nations 


is less than the interaction between their mutual fears, then the arms race 
diverges, 


The Geometry of Stability: An Alternative Approach 


Г the common and both substantively and 
y my “State system, such as that portrayed 
by the classic Richardson model Thi 


t S, and in Chapter 5 we apply it to the 
analysis of a nonlinear System. The i 
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as useful in its own ri Р 
linear КЕ e pm for linear models, but preparatory for the non- 
шске: h solutions in closed form are frequently elusive, and 
ale: ing thus assumes, necessarily. a more prominent analytic 
The conditi ; "T 
Siem кү np eee of this technology is that, except for 
ра ешип ie | oes not occur explicitly in the change equations. 
of інве T be treated as a fixed parameter for purposes 
technique for лаг н. LA This allows us to develop an analytic 
anata attention fa gie e system with time eliminated, and hence with 
slates. Sud à р (where it should be) on the interaction between 
Siewert ст (with explicit dependence on time removed) is 
Tis солы iterature as an autonomous system. 
plane Where bibe ins | is to represent the behavior of the system ina 
Seite the other dod one system state and a second axis repre- 
model the жи гуш state. In the case of our two-nation arms race 
plane: The жю representation, then, is the state space of the (X,Y) 
Sind swe shat pp of system states onabscissa and ordinate isarbitrary, 
iiis omena T j on the ordinate and X on the abscissa. The goal of 
for (1) determirin Й rait in the state space is to provide a graphical means 
ӨЧӨ system fn ШЕ ог locating equilibria and (2) understanding the motion 
Consider th е neighborhood of the equilibrium. 
43. е original arms race model, specified in equ 


ations 4.2 and 


AX, 
AY, 


-aXı + КҮ, +g 


-a'Yı + K'Xı + gl 


form (because no parameter is 


This ; 
System is already in autonomous 
he system is ready for analysis 


directly ti 

y time depend 

ás written. pendent), and hence t 

T ү 
The he system is at equilibrium when th 
We measures of system sta 
d» s through time (AX, = AY: = 0). Thus, at equi 
< May be written as: 


ere is no longer system change: 
tes maintain constant numeric 
librium in X, equation 


КҮ, = aX. - р [4.39] 


and Я 
» Provided k does not equal zero, rewritten as 
Yı = (a/k)X. - g/k [4.40] 
No š А ^ 
ime that Y is expressed as a function of X and fixed parameters with 
єй, аге suppressed. The form of the function is a straight line, and it is an 
т matter to plot this slope-intercept form in the state space plane (X, Y). 
that plot is exhibited in Figure 4.1, giving the locus of points in X and Y 
Produces a steady state in X—8 numeric level in state X that is not 
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(xt Y*) 


Figure 4.1 Phase Plane Representation of Equilibrium in a Two-Nation Arms Race 


3 OT AE rut к ý ritten as 
changing. Similarly, at equilibrium in Y, equation 4.3 may be written А 


Y = (K'Ja)X + gja [4.41] 


This locus of points in the (X,Y) plane is satisfied when state Y is no 
changing in numeric value. 


The geometry of Figure 4.1 is worth careful study. As long as the two 
loci are not wholly coincid 


likelihood 
cepts must be very small. Finally, 


É different slopes produce different values 
for the system equilibrium (X*,Y 
first quadrant, giving a substanti 


Ка < а/к [4.42] 


And, algebraic manipulation of inequality 4.42 Produces an interesting 


result: 


kk’ <a’ [443] 
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Which is identi ; А 
"Thus. тне aia Mas Suc 4.36, the stability condition for the model. 
8 and g' are both s o shows that if the process converges, and if 
equilibrium —an e RE convergence will be toward a first-quadrant 
tation exists. Put ч иш for which a real-world descriptive interpre- 
ducing a first- SHORE ways 1 the geometry in Figure 4.1 holds, pro- 
Conversely, pri ae apt equilibrium, then the equilibrium is stable. 
will not be fried in a runaway arms race is possible in the model, it 
Theapproacht jis with a stable first-quadrant equilibrium. 
by assigning ЖО equilibrium can be studied with this geometric strategy 
using the ériginal conditions (for X, and Y.) апда set of parameter values, 
one-period bod specification of equations 4.2 and 4.3. On this basis a 
expressed in th ment on X and Y is computed, and this movement is 
Жап Y mov € state space as a single vector arising from the separate 
starting point ie eae The tail of the vector is located at the arbitrary 
values), with len Neon of initial conditions not equal to the equilibrium 
assignments ae direction determined by the particular numerical 
and the initial E asan or determined or estimated for the parameters 
values for X and Y попа. 1 his process is repeated, using the computed 
again, as far as o as new initial conditions and then repeated again, and 
the maneria n pleases. The path built up in this fashion describes 
in the (X.Y) st ich the equilibrium 15 approached (or not approached) 
zs ) state space. 
arms "e a theoretical interest in the traditional analysis of 
Ya balance i^ ) a runaway arms race (the future all people fear) and 
situation) Th power (the future we hope for as making the best of a bad 
drawing ч mate cases can be displayed as trajectories ІП the state space by 
illustrate sh € of vectors following the procedure just outlined. We 
Su opeful future of a stable balance of power. ч 
ppose the following assignment of parameter values that satisfy 


Convergence: 
a=.6 a= .8 
к=.9 k'=.4 
g-d g [4.44] 
Then, 
Хе = 1.42 
ўж = 83 [4.45] 
and the equations for the point loci in the (X, Y) plane are 
y = .667X - 11 
[4.46] 


ү = .500Х + -125 
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Substituting the equation 4.44 parameter values into the characteristic 
equation and solving for the roots we obtain 


n = .788 
го = -.428 [4.47] 
and hence the approach to the equilibrium should include a period-two 


oscillatory component. Does our procedure square with these analytic 
results? 


Assume, 
Xo = 1 # X* 
Yo = 1 # Y* [4.48] 


and proceed. The results are plotted in Figure 4.2. The geometry is not 
very interesting but it shows that movement in the state space is toward 
the equilibrium value as predicted. When nonlinear models or less re- 
strictive linear models are analyzed, the trajectories in the state space show 
More interesting patterns, i.e., display more complex motion. This 


becomes evident in the geometric analysis of the nonlinear models in the 
next chapter. 


Feedback Time Delay as a Destabilizing Influence 


In Chapter 2 a stable system is defined asa system generating a sequence 
that approaches equilibrium or returns toward equilibrium after a dis- 
turbance moves it away from equilibrium. An infinite range of models lies 
within the general class of stable systems, producing sequences that 
approach equilibrium in both simple and complex paths, with varying 
degrees of speed. Thus, in an informal sense, some stable systems are more 
“stable” than others, Producing more rapid convergence in a more 
Straightforward fashion. Time delay in feedback loops is an important 
source of these destabilizing influences, and in this section we consider 
its inhibiting effect upon smooth, rapid convergence toward equilibrium. 
Building upon a model by Sprague (1976), the Process of legislative review 
is constructed to examine the consequences of extended feedback 
structures, 

Let us imagine some government agency 
the legislature in order to secure funding. 
within the agency Prepares a proposed agen 
agenda, or demand, serves as a positive inpu 

«as the activity of the agency at time t and 
The legislature, usually an appropriations su 


that makes a yearly trip to 
The bureaucratic leadership 
da of undertakings, and this 
t into the process. (We define 
U: as the proposed agenda.) 
bcommittee, makes a deter- 
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Figure 4. 

e4.2 Phase. Plane Representation of a First-Quadrant Approach Toward 
Equilibrium in a Two-Nation Arms Race 

nda, and agrees to some proportion, 


minati x 
nation regarding the proposed age 
ritten as: 


&, accordingly. This part of the process can be w 

олаш [4.49] 
d with the promise and potential of 
with past performance 


tervals a legislative oversight com- 

ns future agency activity by elimi- 
This review activity, and its 
ve feedback into the system. 


arameter as B, and write the 


Pus. i m is not only concerne 
ànd es ouvi, however, it is also concerned 
miles activity. Thus, at regular int 
Hatin performs a review, and constral 
cubes s. пака practices from the past. 
We ag ent recommendation, serves as negati 
denote this negative review Or feedback p: 
revised model as: 
[4.50] 


Ac = aU, - BAc (к=0,1,2,...) 
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This model is an admittedly crude characterization of the interaction 
between bureaucratic agencies and legislatures, but it serves to illustrate 
several features of the relationship between legislative review and system 
stability. The two parameters are separate and distinct parts of the process: 
a serves as an abstraction of the budgeting process and fj as an abstraction 
of the oversight process. As any student of legislative politics knows, the 
budgeting process includes large doses of oversight activity, but legislative 
reform efforts at creating institutionalized mechanisms for policy oversight 
outside the budgeting process suggest the validity of an analytic separation 
between the two activities. (More realistic models of budgeting are 
developed in the next chapter.) | 

The time subscript for the review process is intentionally left undefined. 
By manipulating k we alter the time delay built into the oversight process. 
When k equals zero oversight is instantaneous, but as k grows larger 
the time correspondence between agency behavior and legislative review 
grows more attenuated. | 

Our goal is to observe the effect of this attentuation. In order to simplify 
our task, we hold о and f, as well as U,, constant while varying k. As 
reasonable values for œ and f, assume that our imaginary agency regu- 
larly wins approval for 75% of its proposed agenda of activities and 
regularly has 25% of its previous activity curtailed through the review 
process. To further simplify our task, think in terms of dollars; We 
measure A, and U, with a dollar metric, assuming that activity reductions 
and approvals proportionally translate into budgetary change. Further 
assume that the cost of the proposed agenda is constant: U, equals 8 million 
dollars for all t. Finally, we trace agency history from three years before 
the instigation of the review procedure, where the three values for A-3, 
A-2, and Ası all equal 6 milion dollars. 

As a first scenario, consider the effect of instantaneous review, k = 0. 
Instantaneous review is produced in very specialized situations. Sprague 
(1976) suggests the example of forest fire appropriations, typically handled 
through supplemental appropriations due to the impossibility of ade- 
quately predicting the number and size of upcoming fires. We represent 


instantaneous review by setting k in equation 4.50 equal to zero and 
substituting values for o and В: 


At = 750, - 25A, = (.75/1.25)0, [4.51] 


Thus, given our constant input of an 8 million dollar budget request, the 
constant output of this process is 4.8 million dollars. More generally. 
instantaneous feedback produces a policy response proportional to 
demand or request, thereby tracking the demand smoothly and without 
delay. 

As a second scenario, assume an annual review so that k = 1. This i5 
typically the situation with a hard-working appropriations subcommittee 
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in the U.S. Congress that performs a thorough review at the same time 
that it considers budget requests. The model becomes: 


A = 150, - 25А [4.52] 


This process, illustrated in Figure 4.3, is oscillatory with period two, 
converging toward a limit that also equals 4.8 million dollars. 
By extending the time delay to two years, equation 4.50 becomes 


A = 750, - 25A: [4.53] 


This form produces complex conjugate roots but is convergent, thereby 
creating a process with sinusoidal-like oscillations of decreasing magni- 
tude” (see part B of Figure 4.3). Finally, by extending the time delay of 
review even farther so that k = 3, we produce the following specification, 


Ay = 750, - .25Ac3 [4.54] 


This time delay results in one negative root and two complex conjugate 
roots, producing a still more complicated time path (see part C of Figure 
4.3). 

Several comparisons among the three specifications shown in Figure 
4.3 are noteworthy. As the time delay increases, the pattern of oscillation 
becomes more complex, and the amount of time required to approach 
equilibrium increases. That is, response time grows longer so that the 
system becomes less efficient. Second, the time delay has no effect on the 
equilibrium value for the system. Thus, by increasing the delay in review 
we increase the complexity of the short-run response, but we do not alter 
the steady-state response. 

What have we learned from 
legislative review can be expecte| 


this exercise? Substantively, delays in 
d to destabilize the policy process. The 
relationship between demand and response is most direct and unambig- 
uous when legislative review is most immediate. If a legislature 15 going 
to exercise oversight, it should not drag its feet. More generally, time delay 
in feedback is destabilizing, if we accept the premise that a process with 
à complex, slow-responding time path is less stable than a process with a 
More direct and faster-responding time path. s 
From a technical perspective such feedback arguments lead directly 
to single-equation models with characteristic equations of higher order. 
The general solution strategies illustrated for the Richardson's model 
are wholly applicable, and the reader might wish to choose a value for k 


in equation 4.50 and attempt a solution. 


Summary 


This chapter has introduced two substantive models: the Richardson 
arms race model and a feedback model of legislative review. These models 
illustrate the frequently occurring necessity of moving beyond first-order 
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Figure 4.3 The Destabilizing Effects of Time Delays in Feedback 


statements of change in depicting social and political processes. The 
solution procedures employed in this chapter are more demanding than 
those employed for the first-order model of Chapter 2, but the solutions 


for second-order models often provide useful analytic mileage. For 
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models of order higher than two, substantively interpreting the solution 
is often quite difficult. As we saw in Chapter 3and will see again in Chapter 
9; some nonlinear models cannot be solved, but a number of alternative 
devices are available to exploit their analytic potential. In particular, 
geometric reasoning plays an especially important role in the two-state, 
nonlinear situation. 


5. NONLINEAR DYNAMIC SYSTEMS: 
BUDGETARY COMPETITION 


This chapter extends previously developed techniques to a consider- 
ation of the budgetary process. The mathematics of this chapter are 
somewhat more demanding because the model that we develop is a system 
of coupled nonlinear equations. Our goal is not to equip readers fully 
with the analytical methods required to conduct similar analyses, but 
rather to illustrate the direction in which the work in this monograph 
naturally moves. 

The appropriations pr: 
American political phenomenon that systematizes and instituti. 
the Lasswellian (1958) formulation of politics as "who gets what, when and 
how," The concern over who gets what in budgeting is the substantive 
problem of this chapter. This question is simply: What (or who) determines 
how the allocation of resources to various budgetary players changes 
Over time? 

The presentation here focuses on agen 
The model that is developed, however, i 
competition over limited resources. The 
ticians, interest groups, fads) competing for a limited set of resources 
(food, votes, favors, members) has quite broad application. The generality 
may be further extended if we allow “predator/prey” interactions in 
addition to pure competition. In such systems, for example, the Lotka- 
Volterra models (May, 1974), one species grows but may be preyed upon 
by other species. Social analogs are obvious. 


ocess is perhaps the most obvious and observable 
onalizes 


cies competing for federal funds. 
llustrates a very general process: 
logic of players (species, poli- 


Incrementalism 


In the area of budgeting one answer to “who gets what" has been termed 
the “theory of budgetary incrementalism." This theoretical approach to 
budgeting derives from studies of public administration and policymaking 
beginning, perhaps, with Simon's analysis of limited rationality in human 
decision making (1957). 

Davis, Dempster, and Wilda 
the incrementalist argument, whic 
Other qualitative research (e... Fenno, 
to their formulation, this year's budge 


vsky (1966) were the first to formalize 
h had been previously supported in 
1966; Wildavsky, 1979). According 
t is determined by the level of last 
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year’s budget plus some increment. Interpretations of this increment 
vary—some suggest it is an inflation factor, others that it is due to the 
inherent nature of a bureaucracy to grow, some that it is merely a descrip- 
tive "as if” theory—but these alternative interpretations are not critical 
to our consideration. The model that seems to best express the incremental 
theory is: 


AAppro, = a(Appro) + error [5.1] 


ог, in more obviously incremental form, appropriations next year for an 
agency are given by all of last year's appropriation plus an augmentation 
or growth percentage. Hence, write symbolically 


Approe = (1 + а)Аррго, + error [5.2] 


The formalization asserts that the change in any particular agency's 
appropriations, in dollars, is some constant Proportion, a (its “fair share”) 
over its current level of funding (its “base” 
tributed error from short-run political 
the process is linear with constant coe, 
of similar budgeting processes fo 
programs. 


). plus some randomly dis- 
and economic forces. Note that 
fficients, first-order, and independent 
r Other agencies, departments, or 


In empirical analyses, œ is usually estimated between .05 and .10. with 
good agreement between observed and theoretically predicted values. 

Indeed, this model and simple variants on it have been utilized to demon- 
Strate budgetary incrementalism at the federal, state, and local levels and 
in the budgetary processes of other political systems (Wildavsky, 1975). 

Notice two consequences of the formalization: (1) There is no equi- 
librium point for the Process, and (2) the dynamics are inherently unstable, 
predicting geometric growth in appropriations over time. These con- 
Sequences may be overlooked if all one desires is a short-term approxi- 
mation of the dynamic Process. The model fails, however, as a descriptively 
adequate statement of the logic underlying the process. 

The incremental models of the budgeting process assume (1) constant 
marginal growth and (2) independence in outcomes for each budgetary 
player. In times when budgetary resources are limited, however, it may 
be more appropriate to treat budgeting asa competitive (not incremental) 
process with interdependent (not independent) outcomes for budgetary 
players. An extended discussion of the empirical justification for such а 
viewpoint and possible methods of model construction is found in Likens 
(1979), 

How can com 


va "i 5 
petition for limited resources be dynamically modeled? 
What are the de 


ductive consequences of such a model? 


Competitive Models of Budgetary Interaction 
From a purely logical point of view, 


à model that expresses change 10 
an agency's or program's share of the bu 


dget should be able to produce at 
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least four distinct qualitative behaviors, dependent on the particular 
substantive political conditions described by different parameters and 
initial conditions. These distinct qualitative behaviors rest on the assump- 
tion that economic conditions are more or less constant. In particular, we 
assume that no radical changes in tax structure occur and that no economic 
miracles or fabulous and unexpected increases in real economic produc- 
tivity are feasible. (But this, of course, depends on the metric of measure- 
ment and some side economic hypotheses. After all, John Stuart Mill's 
steady-state economy is only literary up to this point in history.) 

First, since resources are limited, one qualitative behavior should be 
growth in an agency's share of the budget to some upper limit. The classic 
incrementalist literature takes, as observed data, budget dollar amounts 
measured in current dollars. Taking a share of the budget as a measure of 
agency success provides a sharply different focus. Clearly, raw current 
dollar amounts can grow over time with no change in the delivery of real 
public goods and services if the inflation rate happens to match the 
incremental rate of agency budgetary growth. Equally, even if real growth 
occurs, the growth may be either generalized governmental growth or 
Within budget competitive growth, i.e.. shifts from one agency to another. 
These possibilities can only be considered throughan analysis of budgetary 
shares or shares of some other meaningful economic aggregate such as 
gross national product. Hence, to recapitulate, an agency’s share of the 
budget has some natural upper limit given scarce resources. 

Second, since success is not always guaranteed over time, the model 
should also be able to eliminate altogether some unsuccessful agencies 
Or programs from the budget. Third, since success for some agencies is 
Virtually guaranteed (national defense and Congressmen’ salaries, for 
example, will always secure a share of the budget), the model should be 
able to guarantee a lower limit on some agencies’ shares of the budget total. 
And fourth, it seems likely that many agencies, in the long run, will acquire 
Neither their minimal nor their maximal share of the budget, but will tend 
to attain something in between these two extremes. An acceptable model 
Should, therefore, be able to predict change over time that, in the long run, 
Сопуегреѕ on some equilibrium between the agency’s lower and upper 
bounds. These qualitative behaviors are the most important a priori 
Substantive requirement we impose on the modeling effort. — 

Further, it seems quite likely that the ultimate success or failure of a 
Program or agency will in part depend on its current level of funding. 
That is, a program that is severely pinched for funds will find it difficult to 
Provide services for its (shrinking) clientele; as political support diminishes, 
the chance for budgetary gains diminishes as well (Fenno, 1966). In con- 
trast, an agency that is doing quite will in the budget process can provide 
useful services, expanding its base of political support and maintaining its 
aggressive budgetary posture. This implies, then, that there are multiple 
equilibria toward which an agency's share of the budget will move over 
lime, depending (in part) upon its current share of the budget. 

lt is useful to note that ло linear formalization can adequately represent 
the logic of this process. Since linear systems can at most exhibit a single 
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equilibrium point, a nonlinear representation with the possibility ш 
multiple equilibria is logically required. Regardless of the Sora 
incremental or nonincremental assumptions of the model, the eds we 
properties of the budgetary phenomenon logically require a nonlinea 
dynamic formalization. 


The Logic of Competition 


In order to examine the consequences of competitive interdependence 
among agencies, two nonlinear models are developed in several ied 
First, consider a single agency and denote its budget share by X:. We wis у 
to characterize the changes in the quantity X, in a way that specifies the 


dependence of those changes on the appropriations for a second agency. 
Yı. Hence write 


AX = ПХ, Yi, Ud [5.5] 


Interpret equation 5.3 as asserting that ch 
depend on the agency's past level of succi 
well other agencies 
doing, Ү,; 


anges in the success of agency X 
ess in the competition, Xi; how 
(in particular the principle competitor Y) have been 
and current economic and political conditions that may in- 
fluence relative success rates, denoted by U,. We assume that these 


exogenous inputs, U;, enter the Process in a linearly additive fashion. 
Hence, rewrite equation 5.3 as 


АХ, = ПХ, YJ + g[U.]. [54] 


An analogous argument for a 


gency Y obtains. To keep matters tractable 
we restrict our attention to t 


his bilateral case of two-agency competition 
only. However, it is important to realize that a generalization to multiple- 
agency competitions is Straightforward and follows naturally along the 
lines of the present development. 

The presence of current levels of agency Success, X, and Y, as arguments 
of the function f Tepresents, formally and technically, the substantive 


motivation of our development— agencies are interdependent in the 
budgetary process. 


A number of formulations 
of the system of joint com 
strategy. Define a, 


can be constructed to Obtaina representation 
petition. We follow an operator notation 
and b, as variable rates of change operating on the 
fortunes of agencies X and X respectively. These operators will be assumed 
to function proportionally to the current levels of success for the agencies. 
measured by X, and Y,, and also to be functions of those same current 
levels of Success, X, and Ү,. These assumptions are now sufficient to allow 
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us to write a system of equations for the competition between agencies 
X and Y as 


АХ, = aX, + 0, [5.5] 
AY, = bY, + U, [5.6] 


where the state variables, X, and Yı, denote shares of the budget (total, 
total controllable, total functional, and so forth), acquired by programs 
or agencies X and Y at time t; U, denotes exogenous inputs; а, and Ь are 
variable rates of change for X, and Yı, such that a and b. are both functions 
of X, and Y.. 

Making a, and b, depend on the levels of success X, and Yı, coupled 
With the explicit functional form assumption specified in equations 5.5 
and 5.6, means that the model to be exhibited will be inherently non- 
linear. This is true not because a substantive and formal interdependence 
has been specified between agencies X and Y, but because the formaliza- 
tions of equations 5.5 and 5.6 logically entail nonlinearity given our 
substantive hypotheses on а, and b. namely, that they also depend on 
the current levels of the states X and Y. (In the Richardson arms' race 
model interdependence between states is achieved while simultaneously 
preserving linearity. The present development followsa different strategy.) 

Obviously, it is often the case that more than two agencies or programs 
are in simultaneous competition for resources. Here, however, the 
emphasis will be on the simpler and more easily understood case of only 
two competitors. While some realism may be lost, the practice of studying 
reduced systems is scientifically sound (Hirsch and Smale, 1974: 75). 

Further, to incorporate the central theme of competition between 
agencies and programs, both models include the qualitative assumption 
that an increase in the magnitude of the budget of orie agency decr pio 
the rate of budgetary growth of its competitor (or increases its rate О, 
decay). That is, 


да/дү, <0 [5.7] 


апа 


дь/дх,<0 [5.8] 

Finally, since we are working with budget shares, each Wen" in 
Corporates the assumption that there are levels of funding (in shares 
toward which X; and Y, will naturally tend in the absence chen ee 
and exogenous inputs. Denote these natural limits for X л к оч 
L, respectively, and assume them to be constant for signi ican is 
epochs. This assumption may then generally be expressed as: 


А [5.9] 
4,70 for =, 
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ай. for XL, 15:10] 
a 70 for A. s L [5.11] 
апа 
b >0 fo ү,<1, [5.13] 
T for Yol [5-13] 
b =0 fo Y=L [5.14] 


Expressions 5.9 through 5.14 are abstract and formal, but their ыр 
stantive interpretations should be made clear. For example, p rue e 
5.10 articulates the condition that whenever the current level of bu a 
Support for agency X exceeds its natural limit, the sign of the made 
operator a, should be strictly negative, producing negative rati 
thus moving the state level X, down toward its natural level, Ly, gi = 
that exogenous inputs have been disregarded and budgetary outcomes 
are (temporarily assumed) independent for all players. | n 

There are, of course, numerous formalizations that could incorporate 
these basic assumptions. Two models are studied here. The first р а 
model of “риге” competition in which the intensity of interaction between 
competitors is largely unrestricted. The second is a more realistic repre- 
sentation, where a degree of structural protection exists for each com- 
petitor against encroachment by others. | 

The notion of structural protection is aimed at capturing the idea that 
Some agencies are favored, at least for minimum levels of program effort, 
by the long-run Structure of social, economic, and hence general political 
conditions. Thus, the principal variant of the pure competition model is 


A А : à А e 
one in which a long-run bias Operates in the system, protecting som 
minimal base for each competitor, 


Model I: Unrestrained Budgetary Competition 


To achieve a model of unre: 
acting agencies or programs th 
an explicit assumption for a 
substantively im 
Hence, write 


Strained competition between two inter- 
at captures our central notions, we make 
in equation 5.5. Furthermore, to avoid 
Portant but technically trivial difficulties, we set U: = 0. 


8 = px(Lx - Xi) - сү, [5.15] 


In keeping with our earlier аг 
а, depends on the limit, Lx; 
the current level of success 
duces two new terms in th 


guments, this substitution for the operator 
the current level of success of agency X, Xe and 
of its principal competitor Y, Y,; and it intro- 
€ form of parameters, Ps and cy, all gathered 
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together in an explicit functional form. When equation 5.15 and the 
hypothesis of zero exogenous inputs are applied to equation 5.5 we may 
write 


AX, = [p(L. - X) - 5 Y]X [5.16] 


Equation 5.16 can be read as asserting that change in the level of budgetary 
Success of agency X, AX, is proportional at a characteristic rate px to the 
difference between natural agency limit, Lx, and its current level, X:, 
combined with a discounting factor arising from competition with agency 
Y, operating proportionally to the current level of success, X:. This dis- 
Counting factor operates proportionally to both current agency success, 
X:, and current competitor success, Yı, at a characteristic rate су. The 
minus sign that occurs in equation 5.16 simply indicates the substantive 
notion that competition is at least partially zero-sum--what agency Y 
receives comes partially out of the hide of agency X. An analogous argu- 
ment holds for agency Y. We assume competition is symmetric, and thus 
are justified in writing 
AY, = [p(L, - Yı) - Xi) Ve [5.17] 
To review, the state variables X, and Y: аге proportions ofsome relevant 
budget total received by programs or agencies X and Y. The parameters 
Lx and L, denote the maximal share of the total budget that the agency 
Would receive in the absence of competition, assuming zero exogenous 
inputs. These upper limits (Lx and L,) are assumed to be constant for 
Significant historical periods and are determined by the broad policy 
objectives extant during the era, the general ability of competitors, admin- 
istrators’ political skills, and similar factors. 
The parameters p. and p, denote the genera 
and Y. respectively, as budgetary players. The larg 
is that agency's average net objective performance in 


funding level, Ly or Ly. я age 
Notice that, as the model is written, the greater the agency's averag 


Success, p, the more rapidly it will tend to approach some optimal v 
of the budget. This assumption has been suggested inat least two empirica 
Studies (Sharkansky, 1965, 1968). In this context, an agency ўдава а 
is likely to have greater success as a budgetary competitor if it is heade! 
by People who are highly motivated and articulate political i a 
With experience in bureaucratic infighting. Alternatively, a newly s - 
lished agency or program with relatively unskilled or inexperienced lea: же 
Ship would tend to have a lower level of effectiveness as à budgetary 
Player for some significant historical period. | | d 
Whether ciens motivation, and/or aggressiveness is the pee 
factor is beyond the scope of the model and certainly beyond таан 
available observations. The parameters px and py capture a neto ded 
result of an obviously complex process. Thus, in this instance the necessary 


] objective success of X 
er px or py, the greater 
securing its optimal 
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simplification of reality is especially stark, but without this simplification 
the modeling task could not proceed. Such simplification is justified if it 
yields sufficient insight into the logic underlying the process. ем 

Finally, the parameter c, denotes the rate at which Y competes against 
X, and c, denotes the rate of competition by X against Y. These parameters 
provide a measure of how significantly or successfully each agency en- 
croaches on the other's funds. 

The logic of the model is straightforward. It asserts that, in the absence 
of competition (c, = c, = 0), Xv and Y, approach their upper limits (Ly and 


Ly) according to the logistic law (see Chapter 3) for each player, i.e.. 
independently: 


AX: = p(L. - X)X [5.18] 
AY: = pL, - Ү)Ү, [5.19] 


Assuming, of course, that 0 < Py, Px < І, equations 5.18 and 5.19 pro- 
duce the familiar S-curve typical of diffusion processes (see Chapter 3). 
Note, however, that if Px and p, are unrestricted, the system may be over- 
driven and smooth behavior may disappear. 


For competitive processes, c, and C, are assumed to lie between zero and 


unity, and hence the larger either agency becomes, the greater becomes 
its competitive impact on the other. There is nothing in the structure of 
the process that limits how much the agencies can influence each other. 
In this sense, then, the competition ma 


y be characterized as “pure” or 
"unrestrained." 


EQUILIBRIA IN PURE COMPETITION 


Equilibria for the s 


ystem are obtained as always, by setting AY: = 
AX, = 0: 
= К 20 
O= [p,(L, - X*) - c, Y*] X* [5.20] 
E ус 5.21] 
07 [p (L, - Y*) - c, x*]v« | 


There are four simultaneous solutions for (X*, ҮЗ): 


(X*, Y*) = (0,0) 


[5.22] 
(X*, Y*) - (0, Ly) [5.231 
(X*, Y*) 2 (L.,0) [5.24] 
(X*. Y*) = ([p,p.L. ~eyLyI/[Pyp, = eye], [5.25] 


[p,PyL, = ср; - сус, ]) 
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Net change ceases, then, under four conditions: 


(1) when both agencies are eliminated, 
Q) ^n agency X is eliminated and agency Y obtains its upper limit, 
y. 
(3) when agency Y is eliminated and X achieves its optimal level, 
Ly, and 
(4) when both achieve some competitive level between zero and their 
upper limits. 


Given that three of the four possibilities end in elimination of one or both 
agencies, it is clear that unrestrained budgetary competition is quite 
Darwinian. In addition, this unlimited competition—with its rather 
extreme consequences—is unrealistic in the context of contemporary 
federal budgeting. Agencies and programs in the real world are seldom 
totally eliminated. At worst they tend to move to some minimal level of 
funding that they then maintain year after year. 

This model might well describe budgetary compet 
Settings that are not highly bureaucratized. It may describe, fore 
Programmatic conflict in newly initiated agencies where priorities 
Well established and bureaucratic inertia has not yet mounted. 


ition in political 
xample, 
are not 


PHASE-SPACE REPRESENTATION OF THE NONLINEAR SYSTEM 
Unlike the nonlinear model of Chapter 3, this nonlinear system 
Possesses no closed-form solution. Thus, our choice of analytic strategies 
15 more circumscribed. The qualitative behavior of the model is most 
easily studied in the phase space. That is, rather than examining the be- 
havior of X, and Y, each as a function of time, examine Y: as a function of 
X (or vice versa) through time, i.e., eliminate explicit representation of 
time. (As outlined in Chapter 4 this is the geometric strategy for studying 
autonomous systems.) 
^ A State-space or phase-space representation proje 
Slonal graph of X,, Y;, and t onto the two-dimensional [X:, Y.] plane. Such 
a projection in graphic form is usually termed a "phase portrait The 
Principal question is: Under what substantive conditions are the different 
equilibria reached? What is the likelihood that one or both agencies may 
be eliminated, or that both agencies survive over time? Qualitatively, the 
Phase portrait exhibits a picture of these dynamics that emphasizes the 
Substantive competition of the two agencies sequenced in any arbitrary, 
discrete-time metric. 
"ELS Projections may be displayed as graph 
erent param assumptions. Thus, ; Sta 
the qualitative ае іп т behavior тау Бе easily (and vividly) 


Exhibited in the state space (X,Y). Begin by setting equations = = 
5.17 €qual to zero and obtaining one nonzero point locus for AX; =0 ап 


One for AY, = 0. These nontrivial loci аге 


cts the three-dimen- 


s, which are compared for 
as parameter values vary 


Ү, = [р,/ојх, *[p/o]h (for AX = 0) [526] 
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and 
Yu = [7&/ p]X« + Ly (for AY, = 0) [5.27] 


For this particular nonlinear model it is possible to obtain explicit repre- 
sentations of these zero-change conditions for Y as a function of X. The 
fact that these loci form two straight lines is just well-planned good luck: 
there is no guarantee that more complex nonlinear models will always be 
so obliging. If the model proves intractable at this point, a local stability 
analysis is probably all that is possible. 

For our purposes we disregard the zero loci formed by the X and Y 
axes and focus on those given by equations 5.26 and 5.27. By construction. 
when X; is on its zero-change locus, the change in Xi equals zero; and when 
Y, is on its zero-change locus, the change of Y, equals zero. One implica- 
tion, then, is that whenever and wherever these curves intersect, that 
intersection is an equilibrium point for the system, since АХ, = AY: = 0. 

The question, then, is what happens when X, and Y, move off these 
zero-change loci? Using graph techniques, it is easily determined for this 
model that, whenever X, is to the right of its zero-change line, АХ. is 
negative. We can also determine the behavior of AX, when X; is to the 
left of its zero-change line. Thus, in the first quadrant, at any point to the 
right of the AX, = 0 locus, AX, is negative; to the left, AX, is positive. A 
similar analysis may be performed for the AY, zero-change line, and the 
reader might like to verify that above the AY, = 0 line AY, is negative and 
below it AY, is positive." 

If we then know where (X,, Y.) is with respect to the zero-change lines, 
we can immediately determine the direction of (Xu, Yı) at time (t + 1). In 
effect, each of the zero-change loci "pull" the point (X,, Y.) toward them- 
selves, with Y, able to move only up and down, X, only able to move left 
and right, and with the pull proportional to the magnitude of the dis- 
placement. (A very similar geometric interpretation was originally used 
by Richardson in analyzing his arms race model). Combining the two 
motions gives the vector traveled by (X, Y) from t to (t + 1). For example. 
if the point (X,, Y.) is above ће AY, - 0 locus and to the left of the AX: = 0 
locus, (X:, У) will move down and to the right from t to (t + 1). The reader 
may wish to compare this to the construction in Chapter 4. 


The ДХ, = 0 and AY, = 0 loci produce six distinct geometries in the 
phase space. Three geometries occur when the slope of the AX, = 0 line is 
greater than the slope of the AY, = 0 line. Three different geometries exist 
when the magnitude of the slopes is reversed. The six possibilities are 
exhibited in Figure 5.1, and for each Possibility the qualitative behavior 
of (X, Y.) is displayed. Inspection of Figure 5.1 provides a global analysis 


of the system's dynamics for the first quadrant (the only state space with 
substantive meaning for this model.) 


Algebraic formulations that 


yield direct substantive interpretations can 
be constructed for a number o 


l f these geometries. Recall that p. gives uir 
average success or skill of X and that L, is a measure of the maximal 
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SI :0)Li 
ope for the (AX, *0) Line » the Slope for the (AY, 20] Line 


(вен) 
{ y >(р„/су) 
(py/cy) Ly Ly > (ру/с у Ly <(ру/сх) Ly 


Ly >(p,/cy JL Ly «(o Zell. 


S З 
lope for the (AX,*0) Line «the Slope for the (AY, *0] Line 


Ly <lPx/ey) bx 


Ly <(px/cy)Lx Ly >(Px/cy) bx 
Ly >(Py/ex) Ly 


Lx >(py/ex)Ly L< (py /ex ky 


Figure 5 
Л Geometries of the Zero-Change Loci 


petitors. Thus, the 


s the net objective success of X asa 
L,) measures the agency's budgetary 
es the impact of X on Y, the term 
mpetitive impact of X on the bud- 
tions hold for the terms ру» and 
haracterize various conditions in 
f substantive conclusions: 


Position n 
of X with respect to all other relevant com 


ter: 
wide, may be conceptualized а: 
clout. Si player. In simple terms, (Ps 
гат e the parameter с, express 
gelary 2-4 е interpreted as the net co 
GL. ihera of Y. Similar interpreta 
the first expressions may be used to € 
quadrant, yielding the following set О 
(0 аот X eliminates y if (a) the net objective success of X exceeds 
ak competitive impact of agency Y (pxbx > c,L,). and (b) the net 
jective success of Y is insufficient to defend adequately against 
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pL.» суі, 
cU, S pyly 
рир. «єс, 
У 
Y, 
0 
0 x: 
Figure 5.2. X Eliminates Y 
pb. e, Ly 
Си » p,L, 
Py Px > Cy Cx 
ү, 
% 
X, 


Figure5.3 X Eliminates Y 


competition by X (cL, > PyL,). Typical phase portraits are €x- 
hibited in Figures 5.2 and 3:3. 


(2) Agency Y eliminates X 
That is, the success of Y 
while X lacks sufficient 


8 s ЯГА а. 

if the above inequalities are x 
Overcomes the competitive impact S isl 
clout to defend itself adequately aga 
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pL, < суіу 
ely < руі; 
рур, < Сус; 
Y, 
0 
0 Xt 
XM UL E. 
Figure 5.4 Y Eliminates X 
ba dg EE ect 
p,L. < c, Ly 
cL, € Py Ly 
рурх > Cy Cx 


Figure 5.5 Y Eliminates X 


: = his set of 
encroachment by Y. Phase portraits ШЫ e р 15 
political conditions аге exhibited in dese enough to defend 

(3) Agencies X and Y both survive if both оа, by the other, and if 

ly against encroachment > stabilize the 
See ema on gon ало die ds 
agencies" interaction (PiP: > cc). Figure >- p 
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Pyle > cy Ly 
pyly > AN 
сус, < Py Px 


Xt 


Figure 5.6 Destabilizing Competition and Mutual Coexistence 


pula > суі 
Byly > eL, 
сус, > Py Px 


% 


Xt 


Figure 5.7 Destabilizing Competition and Mutual Coe: 


xistence 


phase portrait of the situation in which the agencies manage mutual 
coexistence. When competition becomes sufficiently intense, how- 
ever, it tends to destabilize even this situation and the result is the 
elimination of one agency or the other. As Figure 5.7 Ls eem 
initial conditions determine the ultimate success, with the relatively 
Stronger agency at time t = 0 finally prevailing. 
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PATTERNS ELIMINATING ONE AGENCY 
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Budget Shares 


Figu == 
ге 5.8 Histories of Appropriations for Competing Agencies as 


of the model that à variety of 
ing agencies. While all of the 
nted here, it is instructive to present 
tterns predicted by the model. 
tic histories that are generated 


e portraits 


"T 
t is apparent from these phas 
he compet 


a be generated for t 

several a series cannot be represe 

Flinies he more common іпбсгаспот ра 

by thi 8 exhibits several typical determinis 
is nonlinear dynamic structure. 
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Model II: Restricted Competition 


Constraints may be imposed on the interaction described by the pre- 
ceding model in a number of ways. We confine our discussion to a single 
elementary modeling strategy—placing lower bounds on X, and Ү,, such 
that spending below these minimal levels is highly improbable. 

Denoting these levels as B, and B,, these lower limits may be incorpo- 
rated in the preceding model by writing 


АХ, = [pL« - X) - GY JXiEX: - By] [5.28] 
AY. = [p(L, - Y) - ХҮҮ, - Bj] [5.29] 


The central logic of competitive interdependence remains unchanged. 
The important difference is that the interaction of X and Y now occurs in 
à structurally limited context. For example, if X is losing ground— moving 
closer to its lower bound—the term (X. - B,) gets smaller and smaller: The 
effect of further interaction is increasingly reduced. And as Y approaches 
its lower bound, By, the impact of X on Y is increasingly minimized. This 
is a modeling device that is frequently useful and worth noting. 

Forcing this constraint on the interaction of X and Y has the technical 
effect of moving the dynamics froma coupled quadratic system to a cubic 
one. The substantive consequences are striking. The most obvious change 
is the multiplicity of outcomes that is now possible. The reader should 
convince himself or herself that not only are the four equilibria of the 
Preceding model possible, but Jive additional equilibria now exist. While 
we do not exhibit these equilibria, their interpretation follows a now 
familiar logic, The System stabilizes when both competitors are eliminated, 
when both are at their lower limits, when One is eliminated and the other 
is either at its minimal level or its maximal level, when one is at its minimal 
level while the other obtains more than its minimum but less than its 


maximum, and when both maintain levels that lie between their lower and 
upper bounds, 


Yı = B, (for AY, = 0) and x, 


always eliminated. 
The model, then, can 


| Benerate all of the qualitative behaviors that are 
logically required of a mo 


del of budgetary competition. Placing limitations 
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on "pure" competitive processes generates dynamics that seem both 
logically and politically realistic. To summarize, the model provides for 
growth in an agency's share of the budget to some upper limit, allows for 
the elimination of unsuccessful agencies from the budget, guarantees a 
lower limit on some agencies’ shares of the budget total, and predicts 
change over time that, in the long run, converges on some equilibrium 
between an agency's lower and upper bounds. 


Conclusion 


The models presented in this chapter attempt to express the dynamics 
of resource allocation as inherently competitive, and often conflictive, 
in nature. It is worth emphasizing that while the analysis does not depend 
critically on the presence of only two competitors, increasing the number 
of players is not without consequence. 

From a technical perspective, increasing the number of states in a 
system, particularly if the system is nonlinear, makes it increasingly diffi- 
cult to obtain a global analysis of the system. With only two states, graph- 
analytic techniques may be quite instructive. With three states, graphics 
are three-dimensional and tricky at best. With four or more states, graph- 
analytic techniques break down altogether, and a local stability analysis 
must suffice. 

From a substantive perspective, increasing the number of competitors 
will produce an increasingly rich political fabric in which a tremendously 
varied dynamic landscape may occur. In the sense of LaPorte (1975) or 
Brunner and Brewer (1971), increasing the number of competitors (and 
hence the level of interdependence in the process) naturally increases the 
probability of generating unanticipated consequences that may or may 
not be socially or politically optimal, and these may be generated by only 
small perturbations of system states, ОГ small alterations in parameter 


values. 


6. DYNAMIC MODELING AS A 
SCIENTIFIC ENTERPRISE 


Dynamic modeling in the social sciences proceeds at many levels of 
abstraction, in widely differing substantive contexts. The basic request 
the investigator makes of any model in any context is this: Does the model 
do the scientific work required of it? Put another way, does the model 
adequately respond to the problem that motivated its development? 
A continuum of possible models ranges from detailed, data-oriented, com- 
plex, descriptive models aimed at practical use to abstract, theoretical, 
non-data-oriented, simplified, only crudely descriptive models aimed at 
providing general insight (May, 1974). We have worked on the side of 
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theoretical models, but they are not purely abstract for we have kept 
empirical interpretations continually in focus. 

The question then emerges: Is this a useful and reasonable strategy? 
To their credit, simple models provide a measure of manageability. For 
example, the mathematical theory for linear forms is complete, equilibria 
are global, and stability may be completely and globally analyzed. More 
complex models are more realistic and capture the details of the process. 
but quickly become unmanageable. These are choices that provide differ- 
ing bargains under differing conditions. It is not better to do theoretical 
modeling rather than more empirically oriented descriptive statistical 
modeling. Nor is it better to do the latter. The models we offer clearly 
possess empirical content with relatively straightforward, substantive 
interpretations. They are, however, far removed from the pole of detailed, 
complex, descriptive models; and from the opposite pole of more abstract, 
theorem-proving, pure mathematical models (Ferejohn and Noll, 1978). 


Thus, the present offering occupies an empirical and theoretical middle 
ground. 


An Overview of Model Structure 


The reader may find it useful to review the logical forms of the models 
we have developed in this monograph, their common substantive inter- 
pretations, and the principal differences among them in terms of quali- 
tative behaviors. The first model of mobilization was a first-order linear 
difference equation with constant coefficients. It may be characterized as 
the time series of a single state variable, in this case voting data, and uses 
only that single observable (and its lag) in model construction and esti- 
mation. Since the model is linear the solution form is known, the equi- 
librium point can be determined, and its stability can beevaluated globally. 


And, even though the model is linear, it potentially exhibits four different 
curvilinear behaviors in the time 


limit, oscillatory convi 
oscillatory divergence. 


domain: monotonic convergence to à 
ergence to a limit, monotonic divergence, and 


The Second model is a first-order nonlinear (quadratic) difference 
equation with constant coefficients, Again, the model uses a time series 
of a single State variable, such as voting data, but in this instance we 
employ the variable, its lag, and the square of its lag. Social interaction is 
built into the model by making a random mixing assumption, and other 
similar nonlinear models have been used to study social interaction 
processes in voting (Kohfeld, 1979) and mass compliance (Likens and 
Kohfeld, 1980; Kohfeld and Likens, 1982). In general, the theory for non- 
linear models is incomplete, solution forms are not usually available, 
equilibria are multiple, and only local stability can be analyzed. The 
quadratic is a special form, however, and results that can be applied to 
determine global stability are available (Chaundy and Phillips, 1936). 
Nonlinear models produce an especially rich variety of time paths: the 
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logistic growth curve, extended period oscillation, and other more com- 
plicated patterns. 

The third model is a system of first-order linear equations that also 
can be written as a single, higher-order linear difference equation. Since 
the higher-order equation is linear, the underlying mathematical theory 
is complete. Systems of first-order linear equations can be used to represent 
complex patterns of interdependence, and a number of strategies are 
available for exploiting their potential: matrix representations of system 
structure, phase space interpretations of a system's qualitative behavior, 
and stability analyses using the system solution. Higher-order models are 
also useful to an analysis of the substantive and technical consequences 
of time delay in feedback mechanisms. 

The final model presented is a system of nonlinear equations. The 
most interesting feature of the nonlinear system construction is the 
structure of reasoning employed to arrive at the final model forms. First, 
the substantive literature (budgeting) was briefly covered. Second, the 
minimum qualitative properties of a reasonable formal model solution 
were developed. Third, a strategy of modeling was employed that, at every 
step, built on these a priori requirements in developing the model form. 
And finally, general mathematical theory was invoked to investigate the 
properties of the model, and to show that the model possessed the desired 
properties; i.e., the model met the problems that motivated its develop- 
ment. In the case of this nonlinear system, the phase-plane interpretation 
of qualitative behavior once again emerges as an important analytic device. 

These models share several important features: structural parameters 
that are constant in time, dynamic properties that logically flow from 
model construction, and, for two of the models, multiple but inter- 
dependent system states. We consider these features before summarizing 


the steps in model development. 


Multiple and Interdependent System States 


feature of the latter two models, and a par- 
f any multiple-state model, is the concept of 
interdependence. Interdependence means that change in one state of a 
social system depends on (and produces) change in other states of the 
system. This is equivalent to the statement that many elements of social 
reality are dynamically interrelated. Change in one social phenomenon 
often has the effect of a pebble tossed into a quiet pond: Itis likely to ripple 
through many other interrelated activities and phenomena. 

While the specific interdependent structures vary across different 
models, the patterns of interaction are always limited to six distinct 
qualitative types. These different patterns of interdependence may be 
represented in a qualitative sign matrix (see May, 1974; Kohfeld, 1976; 
Likens, 1978) with elements that are (+) for positive change resulting from 
interaction, (-) for negative interaction effects, and (0) for no effect. 


An important structural 
ticularly important feature о 
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Effect of State j oni 
+ - ü 
0) 
T id +,+) =) (+, 
Effect ( 
of State - (-,+) (--) (-,0) 
vet 0 (0,4) (0,—) (0.0) 


In biological terms, (0,0) signifies no interaction, (0,5) т yis 
salism, (0,-) is parasitism, (+,+) is symbiosis, EJ is predation, an ian 
competition. Social science analogs are readily identified: The dyn id 
of the Richardson arms race model, for example, is analogous tö a s 
biotic system, and interagency battles for budgetary resources are 
described as a competitive dynamic system. | NE 

Social dynamics are often interdependent dynamics. Change is Es 
isolated; it is interactive. The dynamic interactions that can occur a x 
limited to a very manageable set of qualitative patterns, and all of thes 


А de * ; ivities and 
patterns are potentially observed within various social activities 
processes. 


Structural Parameters 


; i nd 
The models’ structural Parameters control dynamic behavior, а! 
they share several features. Fi 


assumed to be fixed for finite p 
processes is 


pendent duri: 
Second, t 

of the parameters express rates of change for different components of the 

model. These are t 

rates, rates of mob 

so forth. In additi 


those that express 
endpoints. 
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System Properties 


Once an acceptable formalization of the system structure is developed, 
certain of its logical properties become quite important. These properties 
generally concern how and where the process moves across time. More 
precisely, a dynamic formalization of the process leads to questions 
regarding equilibria, stability, and the descriptive qualities of time paths. 

In focusing on equilibria the investigator is concerned with system 
states that remain numerically constant across time. It should be re- 
emphasized that even when the net change of the system states is zero, 
a tremendous amount of change can be occurring within the process. 
People still change their opinions, agencies still compete for dollars, and 
so forth. The numerical values of the system states remain constant because 
all change is balanced, not because change is absent. If the system states 
approach their equilibrium values over time, the process may be charac- 
terized as exhibiting stability. A more rigorous definition is that a stable 
process will return toward equilibrium when its states are perturbed; an 
unstable process will not. 

Combining information about system equilibria and stability, the 
investigator may then develop a strong sense of the system's qualitative 
behavior, the particular history of outputs that a model generates across 
time. Here the focus is on the particular features of the history, revealed 
graphically through a time path of the states of the system for different 
parameters and initial conditions. Alternatively, insights about the 
qualitative behavior of the process may be gained through examination 
of a phase portrait, the simultaneous trajectories of the states across time 
with respect to each other, with time suppressed as an explicit quantity. 

A number of questions emerge regarding the dynamics exhibited by a 
process, Will the process asymptotically approach some limiting value 
over time? Under what conditions will the process oscillate, with what 
periodicity? Can the process reproduce certain patterns—growth and then 
decay to a limit, for example? How, in short, will the states move across 
time, and what central tendencies should be empirically observed? 


Steps in Model Development 


In summary, the two key components in each of the substantive chapters 
are the formalization and its mathematical analysis. An examination of 
the formalizations and their analytics underscores the importance of two 
separate steps in any modeling endeavor. 

The first step begins with the implicit question: “Of what is reality 
composed?" As Quine (1961) informs us, there is universal agreement 
that reality is composed of everything; disagreement is over special cases. 
More precisely, then, the modeling task begins with the construction ofa 
still life portrait—the abstract characterization of a social process at a 
point in time. In each formalization the social process is conceived in 
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such a manner as to reduce a set of complex verbal intuitions to a more 
manageable set of observable (or unobservable) system states. 

Once an acceptable static characterization is in hand, formalization 
proceeds by asking a second question: What are the laws of change? How 
is reality transformed by the logic of the process with which we are con- 
cerned? This second step forces the investigator to make explicit statements 
and assumptions regarding the structure of temporal interdependence. 
The rigor of difference equations, coupled with the analytic potential they 
provide, allows an examination of the consequences that flow from these 
laws of change. In this way we move closer to our goal: an understanding 
of the logic underlying social process and social change. 
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APPENDIX 
Difference Equation Forms 


Difference equations can be algebraically manipulated to produce a 
number of different forms. Three such equation forms are basic and are 
illustrated below. Consider the following difference equation model: 


AY, = ao - ai Yi- [A.1] 


This model is in a recursive form when an isolated left-hand side system 
state is defined in terms of its own lagged values and constants: 
Yea = ao + (1 7 а), [A.2] 
The model is in a canonical form when all system states are on the left-hand 
side, and all constants are on the right-hand side: 
Yor - (1 - an) Yi = ао [A.3] 
Finally, the model is in a homogeneous form if the left-hand side system 
states are set equal to zero: 


Yur - (1) -а)у 20 [A4] 


Matrix Manipulations of the Arms Race Model 


to solve systems of dynamic equations using Cramer's 
rule. Solution by Cramer's rule calls for forming the ratio of two deter- 
minants. The denominator of the ratio is the determinant of the coefficient 
matrix, and the numerator is the determinant of the matrix of coefficients 
after the right-hand side vector has been substituted for the column of 
coefficients whose sequence number from left to right corresponds to 


the system state sequence number from top to bottom. 
Applying this rule, the system in equation 4.11 may besolved as follows: 


It is often useful 


-k 
det ( Е, ' | 
А+а 
ead а. Chy [А.5] 


X,= 
an (ta -k 
Шү ш (Ata) 
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det y x " ) 


ата T NC ma 
e 
ы ( E ene) 


Computing the determinants for X, and Y, we obtain 


_ Ag + ga’ + ke’ [A.7] 
t А? + (a t a') A aa! — kk' 

Е Дв'+ра+ k'g [А.8] 
t 


A? (a t a')A +аа' -kk 


Now treat the left-hand sides as ratios (X, and Y, are implicitly divided 
by one), and cross-multiply obtaining 


Га? + (a + à) A + aa’ - kk]X, = Ag + ga’ + kp’ [А.9] 
[A? + (a + aA + aa - kk]Y, = Ag’ + ра + k'g [A.10] 


The terms involving delta on the right-hand side drop out since the first 
difference of a constant is zero, leaving only well-defined constants whose 


Xun #(2 - a - ауу, (I+ aa^ - kk’ a - aX, = ga’ + кр ГАІ 


Yur + (ata’— 2)Ү + (1 + aa’ - ke — a - a? Y, = g'a + Кв [А.12] 


nical forms used below in obtaining a 
€ system specified in equations 4.9 and 4.10. 
ns А. and A.12 as 


Equations A.11 and A.12 are cano, 
closed-form solution for th 


Finally, rewrite equatio 


Xa (2-a -= a)Xu + (a *a'* kk’ аа 1)Х‹ + ga’ + kg’ [A.13] 
Yo =(2-a - а). + (at are kk’ - aa’ - 1)Y, + gia + k’g [A.14] 
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Equations A.13 and A.14 are identical with equations 4.7 and 4.8 obtained 
earlier by main force and artistry. 

The approach by Cramer's rule is especially convenient for determining 
the equilibrium values for X, and Y,. When the original left-hand side 
deltas are all set to zero, all determinant terms involving delta disappear. 
This leads immediately, using equations A.7 and A.8 to 


X* = (ga' + g'k)/(aa' - kk’) [A.15] 
and | 
Y* = (ва + gk’)/(aa’ - kk’) [А.16] 


as the values of the equilibrium vector S*. 


Continuous versus Discrete Time Models 


The similarities and distinctions between difference and differential 
equation models are best illustrated by considering the solutions and 
behavior patterns of these two types of equations. The analogies between 
the two should not be underemphasized. Existence and uniqueness the- 
orems, theorems on qualitative behaviors, and even the solution forms 
themselves appear in highly analogous ways in the corresponding chapters 
of difference and differential equation texts. These mathematical simi- 
larities, however, should not obscure the fact that there are important 
differences between the two types of equations that are particularly 
relevant in applied work. In some instances solutions can be determined 
for a differential equation, but not for the analogous difference equation. 
Furthermore, even if both equations can be solved, the solutions may be 
quite different: In particular, the qualitative and the quantitative behavior 
of the two equations will be different under certain conditions. To illustrate 
these points we discuss the availability and nature of solutions for first- 
and second-order equations. These equations are chosen because social 
Science applications seldom involve higher-order equations, and the 
extension from a second-order equation to equations of higher-order is 
generally straightforward. " 

The first-order linear differential equation is the simplest to consider, 
and can usually be solved by obtaining an integrating factor. Thus, 
assuming that the relevant integrals can be found in tables of integrals 
and using various techniques of integration, solutions for these equations 
are tractable for most applications. For the first-order linear difference 
equation general solutions depend upon finding inverse difference 
Operators—a difficult task at best. Thus, aside from the special case in 
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which the coefficients are constant, it may be very difficult to solve even 

first-order linear difference equation. А ; 
oem higher-order linear equations and for nonlinear үз Mn 
advantages of differential equations that accrue asa сш a | 
become less significant, although by по means unimportan Е nios 
the second-order linear differential equation can be solved if one a iin 
of the homogeneous portion or a particular solution а сал 
equation is known. Analogous information for the secon {айпа ee 
difference equation renders the problem similar to that 9 and 8. for 
lutions for the first-order equation (see Brand. 1966: chs. 3 an saver: 
details). Thus, neither equation has a readily available soluton. pei 
if information about a portion of the solution is available, it will genera 
be easier to solve the differential equation. "WO 

This conclusion also applies to higher-order linear equations on 
nonlinear equations. In these cases it is often very difficult toifind геа 
tions for either type of equation. There is по general method for so Ss 
even first-order nonlinear differential equations, and the fact that ex E 
tence theorems abound is of dubious practical value. pues ope 
available for special kinds of nonlinear equations, and a TT RE 
techniques can be used to transform equations into different formas Es 
none of this guarantees that a solution for any given equation can x 
found. In nonlinear cases, it is more likely that a differential equation itis 
be solved, simply because more work has been done in this area, but i io 
also true that neither differential nor difference equations can be said 
have readily available solutions. the 

A second distinction between the two types of equations concerns i 
nature of the solutions. Under certain conditions, analogous pairs h 
difference and differential equations will exhibit very different cog 
This can be illustrated by considering a class of equations with wis 
available solutions for both difference and differential equations: use 
equations with constant coefficients. In the case of analogous pairs xai 
first-order linear difference and differential equations, the same En 
conditions and coefficient values can generate qualitatively different pm 
paths. For the difference equation the possible kinds of behavior Ex 
monotonic divergence, monotonic convergence, divergence with osci 
lation, and Convergence with oscillation. For the differential equation. 
the possible kinds of behavior are monotonic div 
convergence. The equilibrium v. 
pairs of difference and differ 
equilibrium point will be the 
derivatives equal zero, Thus, 


ergence and monotonic 
alue will always be the same for analgons 
ential equations since, by definition, t 

point for which successive differences ОГ 


5 A r eater 
the difference equation allows for a grea 
diversity of patterns of qualitative behavior (Samuelson, 1974). 


a . үз HO imilar 
Second-order linear €quations with constant coefficients exhibit p 
patterns. As in the first-order case, two major distinctions occur: (1) 
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ference equations are less stable, and (2) difference equations exhibit 
period two oscillations in addition to the behavior patterns common to 
both equations. This means that the same initial conditions and parameter 
values governing change can result in different qualitative as well as 
quantitative behaviors. 

Difference equations are generally less stable than their differential 
counterparts: An equation that converges in differential form may diverge 
in difference equation form, but the converse is not true. For nonlinear 
equations, conditions for the local stability of difference equations are 
more stringent than the ones for corresponding differential equations— 
the parameter space for difference equation convergence will be a subspace 
of the convergent parameter space of the corresponding differential 
equation (May, 1974: 24-30). Difference equations also exhibit a second 
kind of "instability" in the form of oscillatory motion instead of monotonic 
behavior. Technically, the term "instability" is reserved for divergent 
behavior, but oscillatory behavior may appear unstable in a substantive 
sense when oscillations represent large changes in a relatively short time 
period. Difference equations often exhibit oscillatory behavior patterns 
in regions that, for analogous parameters, yield a smooth pattern from 
the corresponding differential equation. These differences in stability 
conditions, both in the technical sense and in the substantive sense, arise 
as a consequence of the mode of formalization. 

In summary, then, since the solutions of the two versions are different, 
it is likely that numerical time paths generated by analogous pairs of 
difference and differential equations will differ. Furthermore, local 
stability conditions for difference equations are more stringent than those 
for analogous differential equations. This conclusion would also appear 
to be true of global stability conditions, where such conditions exist and 
are known. Finally, for some parameter values, analogous pairs of differ- 
ence and differential equations will generate different qualitative patterns 
of behavior even when stability conditions are satisfied for both equations. 
(For a more complete discussion of the differences between analogous 
pairs of difference and differential equations see Kohfeld and Salert, 1982.) 


NOTES 


1. For treatments of difference equations see Samuelson (1974), Baumol (1970), Gold- 


berg (1958), Cadzow (1973), and Boynton (1980). | x мз, 
2. This distinction between diachronic and synchronic change relies upon à discussion 
contained in Cortes et al. (1974: ch. 1). i A 4 А 
3. Period two oscillation is the most basic oscillatory pattern in which positive an 
negative rates of change alternate in each time period: An increase is followed by a decrease 


is followed by an increase and so on. 
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E 
4. At first inspection these mobilization levels may appear somewhat low. Only AM 
of the population is susceptible to Democratic party recruitment efforts. and the Mee 
converges toward an equilibrium of only 37% Democratic party support. The level iun ke 
zation, however, is defined to the base of all eligibles rather than to all voters. The mean oU; 
County turnout rate for presidential elections from 1956 to 1968 was approximately E tke 
Using this figure as a norm, the pool of potential recruits is 78% of average turnout, ч хә 
mobilization equilibrium is 54% of average turnout. Thus, rather than indicating weakness, 
this analysis points to the strength of the Democratic party in Lake County. fons 
5. We are pursuinga solution strategy that relies on the work of Goldberg ( 1958): of 
useful, compressed outline see Cortes et al. (1974). The strategy developed here is е 
several for solving difference equations. One alternative approach— the z-transform: c H 
forth in Cadzow (1973). Sce the Appendix for illustrations of recursive form, canonical form. 
and homogeneous form difference equations. р К ult a 
6. Readers unfamiliar with complex numbers might find it beneficial to cons 
college algebra text. же 
7. We follow Richardson's convention of assuming that all parameters are Reit 
and several later results depend upon this assumption. It should be mentioned, kowe 5 
that it is at least plausible that some parameters might be negative, and some investigators 
allow for that possibility. first- 
8. A single higher-order linear difference equation can be expressed as a system of bi М 
order linear difference equations. By defining an additional system state (Z.). the reader 
should be persuaded that an alternative form for equation 4.53 is: 


Ac = 7Zi4 + aU, 
Z = ВА 


Indeed, one general formulation of the feedback model isolates feedback as a separate, un- 
observable system state (isolated ina separate auxiliary equation) that is eliminated when the 
System is put into a reduced form (Boynton, 1980: 139). 2 

9. Time delays in feedback are destabilizing in the most general sense as well. Higher- 
order lag structures produce a greater number of roots, thereby increasing the likelihood that 
one of the roots will exceed limits on stability. 

10. The statistical procedures typically utilized in these empirical analyses have been 
criticized by Hibbs (1974) as producing biased, inconsistent estimates, and artificially inflated 
tests of significance. The problem stems from serially correlated error terms in the time series 
as a result of lagging endogenous variables on the right-hand side of the estimating equation. 


ns, especially when written in discrete time, it is 
to break down, since trajectories can "jump 
™ equilibrium li int i to the next. The modeler should, therefore, 
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